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Abstract. We give a brief review of the main results of the theory of elHptic 
hypergeometric functions — a new class of special functions of mathematical 
physics. We prove the most general univariate exact integration formula gen- 
eralizing Euler's beta integral, which is called the elliptic beta integral. An 
elliptic analogue of the Gauss hypergeometric function is constructed together 
with the elliptic hypergeometric equation for it. Biorthogonality relations for 
this function and its particular subcases are described. We list known elliptic 
beta integrals on root systems and consider symmetry transformations for the 
corresponding elliptic hypergeometric functions of the higher order. 
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1. Introduction 

Theory of special functions is widely used in theoretical and mathematical physics 
as a handbook collection of exact mathematical formulae together with the meth- 
ods of their derivation. This concerns the series summation formulae, exactly 
computable integrals, symmetry transformations for functions, differential or other 
equations solvable in terms of "simple" functions, and so on. An impetuous buildup 
of such a database, which was taking place in the XIX century and which was on 
a top of priorities of the mathematics of that time, has changed in the XX century 
by an essential deterioration of interest to special functions, investigations of which 
started to be considered as a pursuit of a secondary importance. Such an attitude 
to investigations in this field was justified by an opinion that all principle types of 
interesting functions "with classical properties" (elliptic, hypergeometric, automor- 
phic, and some other functions) have been found already, and it remains only to 
investigate them in more detail. 

The theory of special functions of hypergeometric type was developing dur- 
ing several centuries, starting from the fundamental results obtained by Euler [1]. 
Gauss, Jacobi, Riemann, Kummer and other prominent mathematicians were con- 
tributing to its foundations. The Gauss hypergeomeric function 2^"i(a, ^; c; a;) is 
a canonical example of functions of such a type. According to the approach by 
Pochhammer and Horn [2], the generalized plain hypergeometric series can be de- 
fined as the sums ^„ c„ for which the ratio Cn+i/cn is a rational function of n. In 
1847 Heine has introduced a g-analogue of the 2-F'i-series 2'Pi{a,b;c;q;x) [3]. The 
general series of such type c„ are characterized by the property that for them 
the ratio c„+i/c„ is a rational function of g", where q is some complex parameter. 
Until recent time only these two classes of functions of hypergeometric type were 
known (including integral representations for them) , and many papers were devoted 
to their investigations. 

Unexpectedly, around ten years ago it became clear that there exist hypergeo- 
metric functions of the third type, which are related to elliptic curves. Such objects 
appeared for the first time within the quantum inverse scattering method developed 
for exactly solvable models of statistical mechanics [4j [5] as elliptic solutions of the 
Yang-Baxter equation [H [7] . As shown by Frenkel and Turaev [8] , these solutions 
(called elliptic 6j-symbols) are expressed in terms of an elliptic generalization of the 
terminating very-well poised balanced g-hypergeometric series loffig with discrete 
values of the parameters. A generalized (l-f l)-dimensional integrable chain similar 
to the discrete time Toda chain was constructed in the paper [9] , and it was shown 
that the same terminating series with arbitrary parameters appears as a particular 
solution of the corresponding Lax pair equations. 

The general formal definition of elliptic hypergeometric series, which was sug- 
gested and investigated in detail in [10], describes them as the series X^n^n, for 
which the ratio c„+i/c„ is equal to an elliptic function of n. Within this scheme, 
the case [8] is characterized by the presence of several interesting structural re- 
strictions upon the coefficients c„. There are certain difficulties in the description 
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of infinite series connected with their convergency. Therefore the general elhptic 
hypergeometric functions are defined by the integral representations [TT] . 

The Meijer function can be considered as the most general plain hypergeometric 
function [12^. It is defined by a contour integral of some ratio of the Euler gamma 
functions. For integral representations of more complicated functions of the hyper- 
geometric type one needs the generalized gamma functions the theory of which was 
developed by Barnes (13j and Jackson [14j more than a century ago. The Jackson 
g-gamma function is necessary for the description of g-hypergeometric functions 
at \q\ < 1. A more complicated function is needed when q lies on the unit circle, 
|g| = 1 [m [m [17l m] . These functions, related to the Barnes gamma function 
of the second order, are actively used in the modern mathematical physics in the 
description of quantum integrable models and representations of quantum algebras 

[IHEniEIlEllEslEl]. 

For the definition of elliptic hypergeometric integrals one needs the elliptic gamma 
function related to the Barnes gamma function of the third order. Importance of 
the elliptic gamma function was stressed by Ruijsenaars [18j . who gave to it this 
name and considered some of its properties. Modular transformations of this func- 
tion were described in |25j. In [llj . a modified elliptic gamma function was built 
which remains also well defined in the case when one of the base parameters lies on 
the unit circle. Other aspects of this function were investigated in [Tfll^ . 

The first exact integration formula, which uses the elliptic gamma function, was 
constructed by the author in 29J . It represents the most general known univariate 
integral generalizing Euler 's beta integral [1]. This elliptic beta integral serves as 
a basis for building general very- well poised elliptic hypergeometric functions. The 
first direction of generalizations consists in increasing the number of free parameters 
entering the integrand, which leads to the elliptic analogues of the functions s+iFs 
[m I30j. In particular, in this way one builds an elliptic analogue of the Gauss 
hypergeometric function and derives its properties |31j . The second direction of 
generalizations increases the number of integrations in such a way that the inte- 
grands acquire symmetries in the integration variables related to the root systems 
[Tl , 32, 33, 34, 35, 36 . One of such generalizations leads to an elliptic analogue of 
the Selberg integral. Over a short period of time, in the papers mentioned above 
and cited below in the list of references, a systematic theory of elliptic hypergeo- 
metric functions of one and many variables has been built. The present review is 
devoted to a brief description of this theory. 



2. Generalized gamma functions 

We use the symbols: Z = 0,±1,±2,...; C - the open complex plane; C* = 
C/{0}; R - the real axis; i = \/— T. 

2.1. The Barnes multiple gamma function. 

The Euler gamma function is a cornerstone of the theory of ordinary hyperge- 
ometric functions [1]. Different analogues of this function with many parameters 
have been considered in the mathematical literature of the beginning of the twenti- 
eth century. The most complete investigations of the generalized gamma functions 
belong to Barnes [13] . As a starting point in his work serves a function generalizing 
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the Hurwitz zeta-function [T] 



1 



c(s, u) = j2 ( 2 i^e(s) > 1, 

(u + nY 

n—O ^ ' 

which reduces to the Riemann zeta-function for u — \. Because of the grown 
interest the Barnes theory was sufficiently widely discussed in the recent literature, 
for instance, in [HI [26l Ell EH [37] . 

Let us take m quasiperiods LOj £ C, which we assume to be linearly independent 
over Z for simplicity (the condition of incommensurability). For s,u £ C the Barnes 
zeta function is defined by the m-fold series 

oc ^ 

Cmis,U;Ul)^ — ft = niLUi + . . . + UmtOrn, 



ni ,...,njn—0 



[u- 



converging for Re(s) > m and under the condition that all cvj lie in one half- 
plane defined by some line passing through the zero point of the coordinates axes. 
Because of the latter requirement, the sequences niwi + . . . + rimUJm do not have 
accumulation points on the finite plane for any rij +oo. It is convenient to 
assume for definiteness that Re(wj) > or lai{ujj) > 0. 

The function Cm(s,M;a;) satisfies the following set of finite difference equations 

Cm{s,U + UJj;uj) - (^,n{s,U;U;) ^ -C,n^l{s,U;uj{j)), j = l,...,TO, (2.1) 

where uj{j) = (wi, . . . , w^-i, Wj+i, . . . ,Wm) and (o{s,u;uj) — u"". It may be contin- 
ued analytically (meromorphically) to the whole complex plane Re(s) < m with the 
simple poles at the points s — 1, 2, . . . , to. The Barnes multiple gamma function is 
defined by the equality 

T„i{u;uj) = cxp{dC,n{s,u;uj)/ds)\^^^. 

It has an infinite product representation of the form 

\ ^ oo , fc 



11' 1 + ^ (2.2) 



Tm{u;uj) V Q 

where ^mk are some constants analogous to the Euler constant. The prime of the 
product means that the point ni = . . . = rim = is skipped in it. In particular, 
the function Ti{u;lu) is directly related to the Euler gamma function T(u), 

ri(u;c.)==-=r - 

which can be checked by straightforward manipulations with the Hurwitz zeta func- 
tion. We note that Barnes used in [13^ a different normalization of the r,„-function, 
in which one has ^„iO = 0. 

The function r,„(u;a;) satisfies to finite difference equations of the first order, 
obtained by differentiation of equalities (|2.1|) at the point s = 0: 



Tmiu + ujj;uj) = j-—^^Tjn{u;uj), j = l,...,TO, (2.3) 



where ro(u; lu) 



1 

„i_i(u;w(j)) 
,-1 
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The following integral representations for the Euler gamma function are well 
known 

I 



^' Jo 2sin(7rs)7c« 
where in the first case Re(s) > 0, and in the second expression | arg(— f)| < tt and 
the Hankel contour Ch starts and finishes near the +00 point, turning around the 
half-axis [0, 00) counterclockwise. One can write with their help 

1 f°° t^^^e~"* 



ns)Jo nr=i(i-e-'^^*) 

ir(l-s) f (-t)^-ie-''* 



2- 7c„nLi(l-e--*)^* 

and analytically continue this function in s to the whole complex plane. Using the 
latter expression, Barnes has derived the following integral representation for the 
multiple gamma functions 

where 7 is the Euler constant. 

The values Cm(0,'u;w) are expressed in terms of the multiple Bernoulli polyno- 
mials 

(-1)" 

Cm(0,U;w) = —Bm,m{u\^^), 

ml 

which are defined by the generating function 



} Bm,n{u\uJi,...,UJm)—r- (2.5) 



We shall need in the following the first three diagonal polynomials 
Bi,i(u|a;i) = ^ - 1 

I \ \ 1 / 2 / ^ Ujf+LOo U!iUJ2\ 

B2,2(w|wi,a;2) = [u - (wi +(J2)u-\ \ — , 

-B3,3(wiwi,w2,w3) = — - — - y^'^fe + ^ [yz^k+^y^'^j^k 



k=l \fe=l j<k 



\fc=l / ]<k / 



Theory of the plain hypergeomctric functions is built with the help of the Euler 
gamma function or Fi (u; wi); the g-hypergcometric functions are tied to r2(u; Wi, ^2), 
and the elliptic hypergeomctric functions "live" at the level of the Barnes multiple 
gamma function of the third order, respectively. 

2.2. The elliptic gamma functions. 

Let 0Ji,(jj2,(jJ3 denote complex parameters linearly independent over Z and lying 
in the right half-plane. We define with their help the base variables p,q,r G C: 

27ri— 27ri— 2iri— /o £y\ 

q = e "2, p = e "2, r = e (2.6) 
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and their modular transformed (t —>■ — 1/t) partners 

— 27r'i— ~ — 27ri— ~ — 27ri— n\ 

q — e ^1, p = e , r — e "3. [^■') 

For IpI, |g| < 1, the infinite products 

00 00 

{z;q)^ = Y{{l- zq^), {z;p,q)oo= Ylil-zp'q'') 

j=0 ],k=0 

are well defined and satisfy g-difference equations 

{z;q)oc 



{qz;q)c 



l-z ' 



{z]q,p)^ {z;q,p)^ 
{qz;q,p)oo ^ ^ , {pz;q,p}oo ^ r • (2.8) 

The shortened theta function 0{z:p) 

e{z;p) {z;pU{pz-^;pU = ^(-1) V^'^-i)/^^'^ (2.9) 

plays a key role in our considerations. It obeys the following simple symmetry 
transformations: 

e{pz;p) - eiz-^;p) = -z-^e{z-p), (2.10) 
and has zeros, 0{z;p) = 0, at z = p'', fc e Z. Evidently, 9{z; 0) = 1 — z. For fc > 0, 
we have 

oip^z;p)^^i^, eip-^z,p) = tlM^. 
{-zYpy^i py 2 ) 

We shall need the r modular transformation rule for 9{z:p), for the 

description of which it is necessary to use the exponential parameterization of vari- 
ables: 

^(e~'"=T;e~'"*^) = e"^^'^("l"i''"=)0 (e'"^ ; e'^'^J^ ^ (2.ii) 

where i?2,2(w|cL;i, C1J2) is the second Bernoulli polynomial. In the following it is 
convenient to use the compact notations 

0{ai, . . .,ak]p) := 0{ai;p) ■ ■■0{ak;p), 9(at^'^\p) := d{at\p)0{ar^\p). 

The simplest gamma function can be defined as a special meromorphic solution 
of the functional equation f{u + uji) = uf{u). Following Jackson's approach 14J, 
we shall be connecting g-gamma functions with the meromorphic solutions of the 
equation 

/(u + c^i) = (l-e2-W-2)j(y)^ (2.12) 
where q = g^^^'^i/'^a Introducing the variable z — e^'^'"/'^^^ equation can be 
replaced by f{qz) = (1 — z)f{z). For |g| < 1 its particular solution, analytical at 
the point 2; = 0, is determined by a simple iteration. This yields the standard q- 
gamma function 7g(z) = 1/(2;; q)oc (which can be considered also as a g-exponential 
function [3]). This expression differs from the Jackson q-gamma function 

r^'H-)-7%%(i-'z)'""' (2-13) 

satisfying the equation T^q\u + l)/r^"^'(u) — {1 — q^) / {1 — q), by a simple change 
of the argument and by a simple multiplier. The limiting transition to the ordinary 
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gamma function has the form hm^^^i Tq''\u) — T{u) [l][3j, but for a simphfication 
of g-hypergeometric formulae it is more convenient to use the function jqiz). 

The modified g-gamma function, which remains well defined at \q\ = 1 as well, 
has the form 

where the contour M + zO passes along the real axis turning over the point a; = 
from above in an infinitesimal way. This function appeared in the number theory 
[15] and in the theory of completely integrable systems [T71 HSl [TH] . It figures in 
the literature under the different names: "the double sine" [TH], "the non-compact 
quantum dilogarithm" [5T] , "the hyperbolic gamma function" [TB] . 

Let Re(a;i), Re(w2) > 0. Then integral (|2.14p is convergent for < Re(u) < 
Re(a;i + 102)- Under appropriate restrictions on u and a;i_2, integral (|2.14p can be 
computed as a convergent sum of residues of the poles in the upper half-plane. For 
Im(wi/a;2) > this leads to the expression 

^("'"""^^^ (2-15) 

which is continued analytically to the whole complex plane of u (this expression 
satisfies equation (|2.12[) in an evident way). This g-gamma function serves as a 
main "brick" in the construction of analytical g-hypergeometric functions at |g| = 1, 
which were not considered in the literature until the recent time. 

The modified g-gamma function is proportional to a ratio of two Barnes gamma 
functions of the second order. The general relation of such a type is derived with 
the help of integral representation (|2.4p and has the form [57] 

where Re(a;fe) > and < Re(M) < Re(X]feLi ^fc)- EU there were obtained also 
infinite product representations of these functions analogous to (|2.15|1 . which we 
are not describing here. 

Already in the works of Barnes it was noticed that the Jacobi 6'i(u|r)-function 
can be decomposed as a product of four multiple gamma functions of the second 
order with different arguments. The exact form of such a relation is (see, for 
instance, [26] ) 

-7TiB2,2(u\l,T) 

r2[u; l,T)r2(l + T -u; l,T)r2(u - t; 1, -T)r2(l - u; 1, -r) 

where p = e^'^*'^. The general relation between multiple gamma functions and 
infinite products of the Jackson type has the form [26] 

°° „-7i-iC™+i(0,w;l,a) 



exp 



ni,...,«„=o r™+i(M; 1, a)r™+i(l - u; 1, -a)' 

where a = (ai, . . . ,am), O = niai -f . . . -1- rimam, and Im(Q;j) > 0. 

Following the logic of definitions of the g-gamma functions, we connect the elliptic 
gamma functions with meromorphic solutions of the finite difference equation 

/(ix + ..i) = 0(e2™/-^p)/(u), (2.17) 
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which passes to ()2.12p for p ^ 0. Using the factorization (12. 9|) and equahties (12. 8|) . 
it is not difficult to see that the ratio 

r(z;p,g) = ^ fr ' - ^-'P'^'f ^ (2.18) 

where \p\, \q\ < 1 and z e C*, satisfies the equations 

^{qz;p,q) ^ 6{z;p)T{z;p,q), r{pz;p,q) = 9{z;q)T{z;p,q). (2.19) 

Thus, the function f{u) — r(e^'^*"/'^^;p, q) defines a solution of equation (|2.17p at 
|(7|, IpI < 1, and it is called the (standard) elliptic gamma function [18]. Because 
non-trivial triply periodic functions do not exist, it can be defined uniquely as the 
meromorphic solution of the system of three equations: 

f{u + L02) - f{u), f{u + ujs) = 0(e2™/-2 . ^) j(^) 

and equation (|2.17p with the normalization of the solution f{J2'k=i ^fc/2) = 1. The 
reflection equation for this generalized gamma function has the form T{z;p, q)T{pq/ z\i 
1. For p = 0, we have T{z]{),q) = ^q{z). 

The modified elliptic gamma function, which is well defined for |g| = 1 as well, 
has the form [TT] 

G(u;u;) =r(e'^'^;p,g)r(re-^"*^;g,r). (2.20) 
It defines the unique solution of three equations: 

J{u + = 0(e2™/-i . r)J{u), f{u + C.3) = e-"^-^("l-)/(u) 

and equation (|2.17p with the normalization of the solution /(X^^^j^ Wfc/2) = 1. 
It is easy to check [55], that the function 

G{u-Lo) = e-'3'^^-'("l"'r(e-2"^;f,p), (2.21) 

where |?^| < 1, satisfies the same three equations and the normalization as the 
function (|2.20p . Therefore these functions coincide, and their equality constitutes 
one of the laws of the modular transformations for the elliptic gamma function 
related to the 5i(3; Z)-group [5S]. From the expression (|2.2ip it follows, that 
G{u]Lij) is a meromorphic function of u for lji/uj2 > 0, i.e. \q\ — 1. 

Because of the antisymmetry condition 53,3(X]fe=i ~ "1^) ^ ^B^_^{u\uj), the 
reflection formula for the G-function has the form G{a;u!)G{b;uj) — 1, a + b = 
J2k=i^k- For \q\ < 1 in the limit p,r ^ (i.e., Im(w3/wi), Im(a;3/a;2) +00) 
expression (j2.20p passes in an evident way to the modified g-gamma function 
7(it; cji, (^2)- The representation (j2.2ip provides an alternative way for the re- 
duction to this function (such a limiting transition was rigorously justified in a 
different way in [H]). As follows from the results of [55|, for a fixed domain of 
values of parameters the function G(u;w) converges in this limit to 7(11; wi, ^2) 
exponentially fast and uniformly on compact subsets of this domain. This result 
is important for a rigorous justification of the corresponding degeneration of the 
elliptic hypergeometric integrals. 

Using the theta function factorization (|2.16p , one can consider equation (|2.17p as 
a composition of four equations for r3(M; w) with different arguments and quasiperi- 
ods. This permits to represent the elliptic gamma function as a ratio of four Barnes 
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gamma functions of the third order [26] 



r3(l + (J + T - u; 1, cr, T)T'i{u - a - t;1, -a, -r) 
For < Im(w) < Im(T + cr), one has the representation 

r(e-»; e--, e^-) = exp f-1 f P^^f^-^] , 

\ 2^-^ A; sm(7rA;r) sm(7rfccr) / 

through which the eUiptic gamma function appeared impUcitly in the work of Baxter 
on the eight- vertex model [6] (see also [5| [25]). 

3. The elliptic beta integral 

As a first example of elliptic hypergeometric functions we describe the elliptic 
beta integral, which was discovered by the author in |29j . 

Theorem 1. We consider six complex parameters tj, j — 1, ... ,6, and two base 
variables p and q satisfying the constraints \p\, \q\, \tj\ < 1 and Jlj=i = TO (^^^ 
balancing condition). Then the following equality is true 

— r(z±^-v o) — ntjtk;p,q), (3.1) 

where T denotes the unit circle with positive orientation and 

{p;p)oo{q;q) oc 
= :r~- ■ 

Here and below we use the compact notation 
r(ai, . . .,ak;p,q) := r{ai;p,q) ■ ■ ■T{ak;p, q), 

r{tz^'^;p,q) := T{tz;p,q)r{tz^'^;p,q), r(2;=^2;p,g) := riz'^;p,q)riz~'^;p,q). 

Proof. We take variables z,q^p £ C, < 1, and five complex parameters 

tm,m = 1, . . . ,5, and compose the function 

' " ■ ■ ■ ' " ^(z±^ q) ni<™<.<5 ntmt.;p, q) ' ^'-'^ 

where A — Y[m=i This function has sequences of poles converging to zero along 
the points 

V = Mp\ A-l<z'^+y+l}™=i,...,5,a,6=0,l,... 

and diverging to infinity along their z — > 1/z reciprocals . Let C denote a 
contour on the complex plane with positive orientation, which separates the sets V 
and V^^ (the existence of such a contour is the only restriction on the parameters 
tm)- For instance, for < 1, \pq\ < \A\, the contour C can coincide with the 
unit circle T. Let us prove now that 

, , dz Ani 
, p(z,ti,...,t5)— = - — — , (3.3) 

c z [q;q)oa(p;p)oo 

where from the needed formula will follow after the substitution A = pq/t^. 
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The first step consists in the derivation of the foUowing g-difference equation for 
the kernel: 

p{z, qti,t2, . . . ,h) - p{z, ii, . . . , ts) = g{q^^z, ii, . . . , ig) - g{z, ii, . . . , tg), (3.4) 
where 

giz, t„ ...,.,) . ,(., t„ . . . , *^)n:;^^^,(,^,^,,),(..,^.;,)7- ^'-'^ 

After the division of equation p.4p by p(z, ii, . . . , fs), it takes the form 

e{hz,hz-^;p) ^ e{At-^-p) 

e{Az,Az-i;p) 9iht,n;p) 



h9{hA-p) ( Z^U^^lOitmZ-\p) Ul=lHtrnZ;p) 



■,P)Ul=2&ititrr^;p) V d{Az-^;p) e{Az;p) 



(3.6) 



Both sides of this equality define elliptic functions of logz (that is they are 
invariant under the transformation z — > pz) with equal sets of poles and their 
residues. For example, 

left-hand\ _ e{tiA,tiA-^;p) -A- e{At-^;p) 



hm eiAz-\p) f ^^"-^-'l) = '^(ti^.tiA ■,p) 
z^A ^ '^'\ side / e{A^:v) 



0{A^;p) ^(ii^™;p) ' 

with the same result for the right-hand side. Therefore the difference of expressions 
in two sides of equality p.6p defines an elliptic function without poles, that is a 
constant. This constant is equal to zero because equation (|3.6p is checked in a 
trivial way for the choice z = ti. 

We integrate now (|3.4p over the variable z G C and obtain 

I{qtiM.---M)-I{tu...M)^ ^ j^9{zM.---M)'^, (3.7) 

where /(ti, . . . ^t^) = p{z^ti^ . . . ,t5)dz/z, and q^^C denotes the contour C di- 
lated with respect to the point 2 = by the factor q~^. Function p.Sp has sequences 
of poles converging to zero along the points z = {tm.q'^p'', A~^q°'p'''^^} and diverging 
to infinity at z = {t^q~''^^P^'' , Aq^'^^^p^''^^} for m = 1, . . . , 5 and a,b — 0,1, . . .. 
For the choice C = T, it is seen that at \tm\ < 1 and \p\ < \A\ there are no poles 
in the annulus 1 < |z| < \q\~^- Therefore we can deform q~^T to T in p.7p and 
obtain zero on the right-hand side. As a result, I(qti, t2, ■ ■ ■ , t^) — I{ti, . . . , ts). 

Requiring \p\, \q\ < \A\, we have by symmetry in p and q that I{pti,t2, . . . ,t^) — 
I{ti, . . . ,t^). Further transformations ti q^^ti and t\ p^^ti can be performed 
only if they do not take parameters outside of the annulus of analyticity of the 
function I{ti, . . . , tg). 

Let us suppose temporarily that p and q are real, p < q, and ^ q'^ for any 
n,k = 0,1, .. .. Impose also the constraint that the arguments of t^, m = 1, . . . ,5, 
and A'^^ differ from each other. Let C denotes now a contour encircling V and two 
cuts ci = C2 = [{pq/ A)p~'^ ,pq/ A] and excluding their z ^ 1/z reciprocals. 

Then we can make transformations ti — > tiq'^ , k — 1,2,..., until the moment when 
tiq^ enters the interval \tip,tip'^], after which we replace ti tip~^; this does not 
take out needed parameters outside of the intervals ci or C2. In this way we obtain 
I{q'p~^ti,t2, ■ = -^(^1, ■ ■■ ,h) for aU j, fc = 0, 1, . . . such that q^p~^ £ 
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Since the set of such points is dense, we come to the conclusion that / does not 
depend on and, by symmetry, on any tm- 

Alternatively, we can use the p-expansion /(ti, . . . , ts) = X^J^o ^n(ti, . . . , t^)p^ 
and check validity of the equalities In{qti, . . . ,t^) = /„(fi, . . . , t^) termwise. The 
coefRcients /„ are analytical in parameters near the points tm = (the constraints 
on the absolute values of parameters from below appear from the requirement of 
convergency of the p-expansion). Therefore, we can simply iterate the dilations 
ti — !■ qti until reaching the limiting point. As a result, /„ and the integral / itself 
do not depend on ti and, consequently, on all tm- 

We conclude thus that the integral / is a constant depending only on p and 
q. In order to find its value, which is given by the right-hand side of ()3.3|) . it 
is sufficient to consider the limit in parameters tii2 1- In this case two poles 
approach the contour of integration, and it is necessary to deform this contour so 
that it crosses over these poles. Then it appears that in the limit tit2 1 only 
the residues of these poles have finite values and the integral itself vanishes. (This 
procedure is described in more detail below.) After proving the integration formula 
in a restricted domain of parameter values, it can be continued analytically to the 
domain permissible by the contour of integration C. □ 

There are many ways to degenerate the elliptic beta integral. In the simplest 
case it is necessary to substitute te = pq/ti . . .t^, use the reflection formula for 
r{z;p, q), and take the limit p ^ 0. After this, the elliptic beta integral degenerates 
to the "trigonometric" g-beta-integral of Rahman [39] (which is connected with the 
integral representation for a 8<P7-series }40) ) : 

(g;g)oo f {z'^;q)ao{z~'^\ q)oo{Az] g)oo(Az~^; g)oc dz 

"^^^ Y{m=l{in,Z;q)oo{tmZ-^;q)oo ^ 

_ X\m=\^-^^m'^q)'x> 

OO 

where A = Y[m=i ^ ^- Further simplification of this equality by taking the 

limit ^5 leads to the famous Askey- Wilson integral 

(g;g)oo f (z^;g)oo(z~^;g)oo dz _ (tlt2^3^4; g)oo 

47ri JTY[m=ii't^nz;q)oo{t,nZ~^;q)oo ^ ni<m<s<4(*™*s; g)oo ' 

which serves as a measure for the Askey- Wilson polynomials [H] - the most general 
orthogonal polynomials obeying the classical properties. The first proof of formula 
was based on an elliptic extension of the Askey approach [42] to computation 
of the Rahman integral. Here we presented the proof obtained in the paper [43] . 
which generalizes the method of computation of the Askey- Wilson integral from 

El- 

If we express in the given g-beta integrals infinite products (a; (7)00 in term of the 
Jackson g-gamma function and pass to the limit q ^ I, then we obtain "rational" 
beta integrals over non-compact contours containing the ordinary Euler gamma 
function [J. Their further simplification by special choices of parameters leads to 
the classical Euler beta-integral 

B{a,f3)= /'\"-i(l-t)/3-irfi=£(^ffi), Re(«),Re(/3)>0. 
Jo r(a + (3) 



12 



V. P. SPIRIDONOV 



The elliptic beta integral ()3.ip gives thus the most general (from known ones) 
exact univariate integration formula including into itself the Euler beta integral as 
a particular case. 

Corollary 2. The following Frenkel and Turaev summation formula [8| is true: 



where titc, = (7^^, nm=o ^"^ ~ 1' '^'^'^ compact notation 

k 

9{tu...,tk)n ■■=l[0{tj)^ 

is used for products of the elliptic Pochhammer symbols 

r{tq";p,q) 



j=o r(i;p,9) 



Proof. We replace in integral ()3.ip T by a contour C separating sequences of the 
poles z — tjq'^p^, j — 1, ... ,6, a,b — 0,1, .. . converging to zero from their z — > 1/z 
reciprocals going to infinity. This permits us to remove the constraints \tj\ < 1 
without changing the right-hand side of (|3.ip . Substitute now t^ = pq/A, A — 
nL=i ^m, and suppose that \tm\ < 1, m = 1, . . . , 4, \pt,i\ < 1 < \tz\, \pq\ < \A\, and 
that the arguments of tm, m = 1, . . . ,5, and p, q are linearly independent over Z. 
Then the following equality is true |32j : 

/dz f dz 

AE{z,t)— = K AE{z,t)— + co(t) ^»(^)' (3-9) 

^ ^ |t5<j"|>l,n>0 

where Ai;(z, i) = nl=i r(im^^^P, g)/r(z±^ ^2±i;p, g) and 

Ut=intmti';p,q) _ ^(^i?'";^) A ^(i™i5)n „ 



^(^5-^A^f ;p,g) ' e{tl-p) JLL^e{qt;n%)^^ ' 

Here we introduced the new parameter t^ with the help of the relation J^^_q — q- 
The multiplier k, is absent in the coefficient cq because of the relation limz_,i(l — 
z)T{z;p,q) = l/{p]p)oo{q]q)oa and due to doubling of the number of residues (the 
latter follows from the symmetry of the kernel z z~^). 

In the limit t^t^ q~^ , N — 0,1, . . ., the integral on the left-hand side of p.9p 
(coinciding with p.ip ) and the multiplier co{t) in front of the sum of residues in the 
right-hand side diverge. But the integral over the unit circle T on the right-hand 
side remains finite. After dividing all the terms by co{t) and passing to the limiting 
equality we obtain the summation formula p.Sp . which was obtained for the first 
time in by a completely different method. □ 

Other proofs of formula are given in |1|151 gSl SZl |H iS] . For p and 
fixed parameters, equality (|3.8p is reduced to the Jackson sum for the terminating 
very- well-poised balanced g-hypergeometric series [I]- The left-hand side of 
formula (|3.8p represents thus an elliptic analogue of this g-series. 

Using the modified elliptic gamma function, it is not difficult to construct the 
modified elliptic beta integral [SS], one of the base variables for which can lie on 
the unit circle, say, \q\ — 1. 
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Theorem 3. Let Im{L0i/u)2) > 0, Im{u!3/LUi) > 0, Im{Lo^/uj2) > 0, and six 
parameters gj £ C, j — 1,...,6, satisfy the restrictions Im^gj/ajs) < and 

nt^o ^ " II G{93+9m;^), 3.10 



l<j<m<6 



where 

{q;q)oo{p;p)oo{r;r)c 



2(<7;9)oo 

Here the integration goes along the cut with the end points — W3/2 and UJ3/2 and 
the convention G{a ±b;uj) = G{a + b; Lu)G{a — b; lo) is used. 

Proof. We substitute relation (|2.21|) into the left-hand side of (|3.10l) and obtain 



6 



Re^'"/^ / ^, . „ — . (3.11) 



where a — 2i?3.3(0|cL') — B3^3{gj\uj). Taken restrictions on the parameters 

permit us to use formula p.ip with the substitutions 

z ^ J e % p — > e "^3 ^ q ^ e "3 , 

which yields for p. lip 

2KW3W^^e'^"'/3 



. 9i +9n 



{f;rUip;pU i</<i<6 



(r;r)^(p;p)^ i<,<™<6 



n G{gj + gm;uj), 



where 6 — J2i<j<m<6 ^3,3(5j + 3m|^)- A straightforward computation shows that 

a + b 

Therefore for the choice 



2ti;c,eS(ELi '^'=)(ELi "fc ^) 
j^-^ — 

^^2if;f)ooip;p)oo 

we obtain the needed result. After application of the modular transformation law 
for the Dedekind function 

V / 00 

to infinite products entering definition of we obtain 



iuj2 {r; r)oo{p;p) 

One more application of the relation (j3.12p permits us to replace the exponential 
function by a ratio of infinite products, and this leads to the needed form of k. □ 
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If we take the limit lm(a;3) — > cx3 in such a way that p,r 0, then the modified 
elUptic beta integral reduces to a g-beta integral of the MeUin-Barnes type. More 
precisely, for aji^2 such that Im{uji/u!2) > and Re(wi/a;2) > 0, we substitute 
ge = X)fe=i ^k—A, where A = 9j ^^'^ apply the inversion formula for G{u; uj). 

After that we set uj^ = ituJ2, t +oo, and obtain formally 

r rij^i 7(gj ±u;uj) du ^ ^ (g; g)^ ni<j<m<5 7^ + gm; t^) 
Jl "fi±2u,A±u;Lu) LU2 (g;<7)oo gj;uj) 

where 7(u; wi, 1^2) denotes the modified g-gamma function, and the integration is 
taken along the line L = iuj2M.. This resuh is true provided the parameters satisfy 
the constraints Ke{gj/uj2) > and Re((^ — uji)/uj2) < 1. This integration formula 
represents a "hyperbohc" analogue of the Rahman integral; it was proved for the 
first time by Stokman in [50j . Because of the non-compactness of the integration 
contour, the described method of derivation of (|3.13p is rigorous under the condition 
of uniform convergence of the function G{u;uj) to ^(u]u;i,lu2), which follows from 
the results obtained by Rains in [28]. One can establish also formula (|3.13[) by 
the method, which was used above for proving the elliptic beta integral. The limit 
q ^ 1 leads to the same rational beta integral as the "trigonometric" Rahman 
integral. 

Summarizing the consideration of the present section, we see that the elliptic 
beta integral includes into itself the whole hierarchy of exactly computable integrals: 
two types of the g-beta integrals, the rational class of beta integrals, whose kernels 
are expressed in terms of the Euler gamma function, and the classical Euler beta 
integral. This scheme refiects the general picture of degenerations of the elliptic 
hypergeometric functions which was rigorously considered in [28j . 



4. General elliptic hypergeometric series and integrals 

In the papers [TD] and [TT] , the author has proposed definitions of general elliptic 
hypergeometric series and integrals, which will be considered in this section. 



4.1. An elliptic analogue of the Meijer function. 

Univariate contour integrals A{u)du are called the elliptic hypergeometric 
integrals, if the meromorphic function A(u) satisfies the following system of three 
equations 

A{u + LUk)^hk{u)A{u), fc=l,2,3, (4.1) 

where wi.2,3 G C are linearly independent over Z parameters, and hk{u) are some 
elliptic functions with the periods Wfc, uJk+i {we set uJk+3 = Wfe). 

The general elliptic function of the order s with the periods uj2 and can be 
represented in the form (see the Appendix) 

fei(^) = yijy[^ g(^^.^2w^,.^) , 

where yi is an arbitrary constant and tj , Wj are some parameters satisfying the 
balancing condition Jlj^i — Y['j=i^j remind that p = ^'^^^il^'iY Using 
properties of the function r(z;p, g), it is not difficult to build the general solution 
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of the fc = 1 equation in (14. 1|) for \q\ < 1: 

u/cJl / N TT ^{tjZ',p, q) 



^{u) = yT^^{u)\{ 



e 



27Tiu/LtJ2 



where ip{u + lui) = '^{u) is an arbitrary periodic function. So, if we would restrict 
ourselves to a single equation for A(m), then our definition of the integrals would 
be highly non-unique. 

An arbitrary elliptic function of the order i with the periods loi and W3 has the 
form 

h2{u) =2/2 II 



where |r| < 1, 2/2 is an arbitrary constant, and the parameters satisfy the constraint 
n^=i = n^=i "^0- The k = 2 equation from ()4.1|) serves now as a constraint for 
the function </?(it). It is easy to prove that the common solution of these two 
equations for j^/j < 1 has the form 

A(«) = n ^^-Cw^^^'^l n ^^^Tw'^'^l (4.2) 

where 4>{u) is an arbitrary function satisfying the equations (j){u + uji) = yi(j){u) and 
(j){u + UJ2) = y24'{u). 4>{u) is thus a meromorphic theta function with the special 
quasiperiodicity multipliers the general form of which is easily established (see the 
Appendix) : 



^cu+d TT c^lQfce ' -,q) 



fe=l 



where m is an arbitrary integer, the parameter d is arbitrary, and the parame- 
ters ak,bk,c are connected with yi and j/2 by the relations 2/2 = e'^'^^ and yi — 

Due to the representation in terms of the elliptic gamma functions, the function 
(t){u) can be reduced to the pure exponential factor by the replacement of s in 
A(u) by s -|- 2m and the choice of parameters ts+fc = po.kits+m+k = bk,Ws+k = 
ak,Ws+m+k — pbk, k — 1, . . . ,m, which does not violate the balancing condition 
Ojii™ — rijii™ ^j- Since s was arbitrary from the very beginning, we can set 
— e™^'' without loss of generality. 

As a result, two equations determine already the kernel A(u), i.e. for the elliptic 
function h^iu) we automatically obtain 

h.(u) = e'''^' TT ^ ^ ^ TT ^ ^ 

^^3^ 6'(Wje2""/'^2;g) 11 0(r-lf,.e-27r™M;g) ■ 
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To summarize, for |g| < 1 the most general elliptic hypergeometric integral has 
the form [HI E] 



e 



with two balancing conditions for the parameters indicated above and some inte- 
gration contour C . 

Consider the definition of integrals for |g| = 1. It appears that now even the 
function h2(u) cannot be arbitrary. In this case it is necessary to take i — s and 
choose the parameters tj , tj , and Wj , Wj in such a way that all F-functions are 
combined to the modified elliptic gamma functions G{u]uj) (it is in this way that 
this function was built in [H]). This leads to the integrals of the form 

-^^ng^^rfu, (4.4) 



e 



where the parameters gj and vj satisfy the balancing condition X]j=i(5j ~ "^j) ~ ^ 
together with the relations tj = g^'^^fj/'^^ wj = ^ = ^e^^Trig. /^i^ ^ 

re-2™^/'^S and yi,2 = e'^^-^ 

The case I9I > 1 appears to be equivalent to the case |g| < 1 after a change of 
parameters and leads to the integrals 

gcu+rf TT i (,g Wje ) -pr Wje > ,q ,r) 

c j-J-^ r(g-iije2""/'^2;p^g-i) 11 r(g-iije-2™/c^i;g-i^r.) ■ 

Functions (|4.3|) . (|4.4p . and (|4.5p can be called as elliptic analogues of the Meijer 
function, because for some particular choice of parameters and of the integration 
contour C they degenerate to that function 12J. During this degeneration pro- 
cedure, at the intermediate steps there appear various g-analogues of the Meijer 
function, including the cases considered in [51 1. The more general theta hyperge- 
ometric analogues of the Meijer function, for which the kernels A{u) satisfy the 
system of equations (j4.ip with (u) given by arbitrary meromorphic theta func- 
tions, are built in llj; we do not consider them here. 

4.2. Well poised and very-well poised integrals. 

We consider integrals (|4.3p with £ ^ c = d = and replace the integration 
variable by 2: = g27riM/(^2 ^ Until now we did not fix the integration contour C. Let 
us choose it as the unit circle T oriented counterclockwise. As a result, we obtain 
integrals of the form 

dz ^ A ^itkz;p,q) 



/a(z)^, A(z) = n 



^ r{'Wkz;p, q) 

with Ilfe=i = 11^=1 '^k- In the case when the conditions w^tk = pq, k — 1, 
are satisfied, the integrals take the form 

r dr * 

/ A(^)(z)-, A(^Hz) = l[r{tkZ,tk/z;p,q), 



k=l 



and are called well poised integrals. The balancing condition for them takes the 
form 11^=1 = (w)* or I\k=i ^k = fJ-ipq)^^^ with the ambiguity in the sign choice 
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II = ±1. The reflection formula T{a,b;p,q) = 1, ab ~ pq, shows that the choice 
of parameters tjtk = pq plays an essential role, since it reduces the number of 
parameters in A(*)(z). In particular, for = pq the variable tk drops out of 
formulae completely. The function 

is evidently p-elliptic, that is h^P\pz) = h^P\z). Denoting Uk ■= tkZ, Vk ■= 
pqz/tk, A := pqz"^, we can rewrite h^P\z) as 

with the conditions of well poisedness UkVk = A, k = l,...,s, and balancing 

Let us consider all possible p-shifts of the parameters ui, . . . , and A: 

uk^p'^'-uk, X^p^X, nk,NeZ, 

and require that /i^*"^ is invariant under the maximally possible group of these trans- 
formations. The balancing condition leads to the constraint 5^^^^ nk — sN/2. For 
= it is easy to check that h^P'>{. . . ,pua, ■ ■ ■ ,p~^ui,, . . . ; A) = h'^P^ui, . . . , Ug; A), 
i.e. is an elliptic function of all its parameters. The transformations with 

N ^ are more complicated and depend on the parity of the variable s. For odd 
s the integer N must be even. The full symmetry group is generated then by the 
transformations with N = and, say, ni = s, Uk = 0,k 1, and N — 2 yielding 

h'-P'^{p'ui,U2,...,Us;p'^X) ^ A" ^ ^ 
/iW A) "nLi«fe ' 

that is the value of fj, is not fixed. 

As to the even values s = 2m, in this case there are no constraints on N, and it is 
sufficient to consider the transformation corresponding to the choice ni = m, N = 
1, with all other n/j = 0. Then we have 

h^P\p"'ui,U2,...,U2m;pX) ^ A"^ 

h(P^Ui,...,U2m;X) UlZ^k 

Requiring this transformation to be a symmetry, we fix uniquely the balancing 
condition, /x = 1. 

It is not difficult to check that for any = ±1 the following equality is true 

for all i,j e Z. Vice versa, from this condition one can derive the balancing 
condition with /i = ±1. Passing to the limits z ±p~*/^(7~^/^ and using the 
symmetry A('')(z) = A('')(z-1), we obtain A^") (±p*/2^J/2^2 ^ A('>){±pi/^q-i/'^f. 
The straightforward computation yields 

aW(±p'/2^^/2) = AW(±p'/2g-^/2). 
For even values s = 2m and /U = 1, we obtain 

A(2™)(±pV2^j72) ^ A(2™)(±p'/2g-^/2), (4.7) 
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provided the functions on both sides are well defined. The latter requirement is 
satisfied provided we do not hit the poles, that is U ^ ip^/^g''/^, a,b gZ, for all r. 

We call the integrals A{z)dz/z very-well poised, if their integration kernels 
have the form 



i-[z ,p,q) 

This kernel can be obtained from A^'^"^~^^'^\z) by restriction of the parameters 
i2m+7, • • • , i2m+i4 = (±(pg)^/2, ip^^^g, ±pg^/^ ±pq); 

this follows from the reflection formula and the argument duplication formula for 
the function r{z;p,q): 

T{z^;p,q) = T{±z,±q'/^z,±p'/h,±{pqy/^z;p,q). 
The balancing condition takes now the form 

2m+6 

n^fe=MM)"+\ (4.8) 
fc=i 

where we count the sign choice /i = 1 as canonical, since it leads to additional 
symmetries. The kernel aI?^^^\z) can be obtained also from A(^'"+^^)(z) by im- 
posing the constraints i2m+7, • • • , i2m+i2 = {±p^^^q, ipq^^"^, ±pq), since the choice 
t| = pq simply removes the corresponding F-factors. However, we shall not use such 
a reduction, because it changes the sign in the balancing condition: Hfe™!^^ = 
-fl{pqr+\ 

In the following we shall be studying the very- well poised elliptic hypergeometric 
integrals of the form 

with the "correct" balancing condition Hfe^i^^^fc = {pq)"^~^^- These integrals rep- 
resent elliptic analogues of the plain hypergeometric functions rn+iFm- In partic- 
ular, for m = we obtain the elliptic beta integral. The constraints t2m+5 — 
{pqY^^, t2m+6 = —{pq)^^^ reduce these integrals to . . . ,t2m+4) with 

the balancing condition 0^=1^*'^'= ~ — (w)™- Appearance of the "-" sign on 
the right-hand side simply indicates that these integrals should be considered as 
some generalizations of ^-i-i^rn-functions, and not of mFm-i- For example, such a 
choice in the elliptic beta integral yields after taking into account of the relation 

^{-pq;p,q) = ^{-p;p)oo{-q;q)oo the following: 
UUintkz^';p,q)dz 



-I 

Im Jc 



2m J c T{z^'^;p,q) z 

4 

= 4(p';p')oo(?^;5^)oo n r(i,tfc;p,g) JjF(p94;p2,g2)^ 

l<3<k<A k=l 

where nfc=i ^fe = ^^"^ ^he contour C separates sequences of poles converging 
to zero from those going to infinity. For p ^ we obtain a special case of the 
Askey- Wilson integral. 
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The multiplier appearing in integrals' kernel from the very- well poisedness con- 
straints for parameters leads to an interesting property [34j 

AZT'^p'^^'^') - ~f^''AtT'Hp'^^~'^% (4.9) 

that is to the "-" sign in the right-hand side for the canonical choice fi ~ I, which 
is sharply distinct from (14. 7|) . Indeed, for i = or j = we have zeros in both 
parts of the equality, and in other cases we obtain 

iim,_p./2,./2rp;i(z2,2-2) _ e{p-\-i-q)e{p\^-p) 



lim,^p./.,-,/2 rpTg\z2, z-2) e{p--qO-q)e{p^q-i;p) 

which together with the described properties of A*^2™+^^(z) yields the presented 
formula. For i,j^O there appears a non-commutativity of two limits: 

^lm^_,pi/2qj/2 Llvwp (Z) _ y 



,. »(2m+6). N ^ 

lim^_^pi/2g-3/2 IXvwp (Z) 

M2m+14)( i/2 j/2 



although both of them are well defined. The reason for such a contradiction con- 
sists in the use in the very-well poisedness condition of the "forbidden" values of 
parameters leading to poles and zeros of ISP''^^^'^\z) for z = p^l'^q^^'^ . 

As shown in [34], the product ni<j</£<m(^j^fe;P' '?)oo-^^™''(^) ^ holomorphic 
function of parameters € C* . Additionally, the "-" sign on the right-hand side of 
(|4.9p guarantees that the integral /(™)(t) is holomorphic in the points = p~°'q~^ , 
a,b — 0,1,.... If some of the parameters take the forbidden values ±p°/2gfc/2^ 
a, 6 G Z, then the latter property disappears. 

4.3. Series. 

According to the general definition [10], formal series X^nez called the 
elliptic hypergeometric series, if the ratio of neighbouring coefficients c„+i/c„ is an 
elliptic function of n G C. This definition lies in the stream of ideas of Pochhammer 
and Horn which are used for building the plain and g- hypergeometric series [5]. 

As we saw already on the example of the Frenkel-Turaev sum, the elliptic hy- 
pergeometric series appear as sums of residues of certain sequences of poles of the 
elliptic hypergeometric integrals' kernels. Indeed, let us consider the poles of the 
integrand in (|4.3p located at the points u = a + luiu, n = 0, 1, . . ., for some pa- 
rameter a, and denote residues of these poles as c„. For u a + ujin, we have 
A(u) Cn/{u — a — Wirt) -|- 0(1). Now it is not difficult to notice that the ratio 

lim — — — hm iLiyu) — hi(a + ujin) 

u^a+uin A\u) C,i M^a+win 

is an elliptic function of n with the periods 102! and wa/oji, which demonstrates 
the general connection between the integrals and series. 

An arbitrary elliptic function h(ri) of the order s 4- 1 with the periods 102/ loi and 
wa/wi has the form 
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where n£i = nfe=i''^fe- We define the elliptic Pochhammer symbol 9{t)n as 
the solution of the recurrence relation c„+i = 9(tq"';p) Cn with the initial condition 
Co = 1: 

n;Jo^(V;p), n>0 

Then it is easy to deduce the explicit form of the formal bilateral elliptic hyperge- 
ometric series 

with the normalization of the zeroth coefficient cq ~ 1. Choosing Ws+i = q and 
replacing ts+i to, we obtain the unilateral series 

fto,ti,...,ts \ 9{to,ti, . . . ,ts)n n tA^^\ 

\wi,...,Ws J ;^^0{q,Wi,...,Ws)n 

For fixed tj and Wj, in the limit p we have 9(t)n (i; <?)n = (1 ^ 0(1 ^ 
qt) . . .{1 — and the function s+iEg reduces to the q-hypergeometric series 

m 

ta,ti,...,ts_ ^ A _ \^ i'to;q)n---{ts;q) 



+i'Ps\ ■,q;y 



E \bo,q)„ . . . \i,s,q)n „ 



Wi,...,Ws' ' / f-'{q]q)n{wi;q)n---{ws;q)„ 



71=0 

with the condition 11^=0 ~ ? 11^=1 ^j- Parameterizing tj = q"^ and Wj 
the limit q ^ I we obtain the series 



^ /uo,ui,...,Us. \ _ (^o)n ■ ■ ■ „ 
V vi,...,Vs ' J ^n!(i;i)„...(w,)„ 

where (a)„ = a{a + 1) . . . (a + n — 1) denotes the standard Pochhammer symbol and 
uo + ■ ■ ■ + Us — 1 + vi + . . . + Vs. The latter constraint for the parameters is not 
essential, since it disappears already for the sf^-i-function obtained after taking 
the limit Us oo. 

Investigation of the conditions of convergence of the infinite series (j4.1ip repre- 
sents a serious problem and it was not solved completely to the present moment. 
Therefore in the applications of the i?-series it is usually assumed that they termi- 
nate because of the condition tk — q^^p'^^ for some k, where iV = 0, 1, . . . , M E Z. 
It is worth of noting that the formal 2-E'i-series does not represent a natural ellip- 
tic generalization of the Gauss hypergeometric function, because it does not obey 
natural analogues of many important properties of the 2-F'i-function. 

Series ()4.1ip are called well poised, if the following constraints on the parameters 
are satisfied taq = tiwi = . . . = tsWs- The balancing condition for them takes the 
form ti- ■ - ts ^ ztq^^^^^/'^tQ^ and the functions h{n) and s+iEs become invari- 

ant with respect to the transformations tj — > ptj, j = l,...,s — 1, and to p^to 
(for this it is necessary to count ts as a dependent parameter). In the same way 
as in the case of integrals, for odd s and the sign choice in the balancing 
condition there appears an additional symmetry — the functions h{n) and s+iEg 
become invariant with respect to the transformation tg —f pto (with the compen- 
sating transformation ts — > p^^'^^-'/^is)- The function h{n) becomes thus an elliptic 
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function of all free parameters logtj, j — 0, . . . , s — 1, with equal periods, that is 
there appears some kind of "total ellipticity" [lOl [31] . 

The next structural restriction, which is needed for building the most interesting 
examples of the series, looks as follows 

ts-3 = qVto, ts-2 = -qVh, ts-i = q\/to/p, ts = -qVpto 

and is called the very-well poisedness condition (it is related to doubling of the 
argument for theta functions). In view of the importance of the very- well poised 
elliptic hypergeometric series, there is a special notation for them [52j : 



— s+lEs 




ta,ti, ■ ■ ■ , ts^i, q^/to, — g-^ 
,qto/ti, . . . , qto/ts-i, Vto, ■ 

where Hfe^i = ±4'' '^^^^9^'*^^^^^ (in this balancing condition for odd s the 
sign choice is considered as canonical). The choice tj = ±y/qto removes completely 
this parameter and reduces the indices of the series by 1 (this choice can change 
the sign in the balancing condition). For sufficiently large values of s, one can 
choose four parameter values in such a way that the very-well poised part in the 
series coefficients cancels out, and it becomes again only the well poised series. For 
y = 1, this argument is dropped in the notation of s+iT^-series. In this scheme, the 
Frenkel-Turaev formula p.Sp yields a closed form expression for the terminating 
ioVg(to;ii, ■ • ■ ,h;q,p)-seTies. 

In order to consider modular transformations it is necessary to use the param- 
eterization tfe — g"*", Wk — q'"'' , q — e^^^" , p = e^^^^ and replace 0(g°;p)-functions 
by the "elliptic numbers" [a] = 6*1 (era |t). Suppose that the parameter y does not 
depend on r. Then, it is not difficult to verify with the help of formula (|A.3[) that 
the functions h(n) and, respectively, s+iEs will be modular invariant under the 
restriction X]fc=o '"fc ^ ^ + Y%=i '^t defined modulo iTja"^ . For completeness of the 
description, we present an explicit expression for the most interesting s+iVs-series 
with odd s in this notation: 



oo r„ 1 2m~5 r 1 

2mV2m-i[q j9 ,q,p,y) — / ^ — f — ] — 11 77—, T-y 

fo ["O] fc=0 [^+^0-Uk]n 



m 



where [a]„ — [a][a + 1] . . . [a„ — 1] and the balancing condition X]/c=i 
4 -|- (m — 3)wo- This series is automatically modular invariant. In [5], the symbols 
s+iCs and s-i-i'^s were suggested for the additive system of notation for series, but 
we do not use them here. 

Importance of the balancing condition for the plain and g-hypergeometric func- 
tions was known for a long time, because in its presence there appear some addi- 
tional identities [3]- The same is true for the notions of well poisedness and very- 
well poisedness. However, the corresponding constraints on the parameters were 
appearing in a rather ad hoc manner, and their deep meaning was missing. The 
elliptic hypergeometric functions clarify the origin of these old concepts. Namely, 
the balancing condition is connected with the condition of double periodicity of the 
main elliptic function used in the construction of series or integrals. The condition 
of well poisedness is connected with the condition of ellipticity in all parameters 
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determining the divisor of this eUiptic function. The condition of very-weU poised- 
ness is tied to the rule of doubling of the argument of theta functions. Strictly 
speaking all these notions are defined in fact in a self-contained manner only at the 
elliptic level. Indeed, there are limiting transitions from the elliptic hypergeometric 
identities involving s+iT4-series to the g-hypergeometric relations such that there 
appear basic s+iyjg-functions which do not obey one of the mentioned properties 
[SSI [52]. 

The degeneration of s+ii^s-series described above at p ^ leads to s+if^s-series 
with the constraint for parameters resembling the old g-balancing condition [2 [3] , 
but not coinciding with it. The limit p — > for s+iVg-series with fixed parameters 
leads to the very-well poised balanced s_i<y9s_2-series having their own notation [3]: 

lim s+iVsito;ti, . . . ,ts-i] q,p; z) = V " — TT " (gz)" 

p^o I -to I^J-^(qto/tk;q)n 

=: s-iWs-2{to;ti, . . . ,ts-i]q,qz). 

In a remarkable way, the balancing condition for the s+iKi-series coincides in this 
case with the usual balancing condition for the s-iW^s-2-series [3j[T0]. The possi- 
bility of fixing the sign in the balancing condition for odd s from the requirement 
of existence of an additional symmetry strengthens the "elliptic" point of view on 
the functions of hypergeometric type and indicates on an indispensable connection 
of these two classes of functions. 

In conclusion of this section, let us mention that the quadratic transformations 
for the s+iV^s -series and related summation formulae were considered in [49l [521 153 1 
I54j . Other specific elliptic hypergeometric series were investigated in the papers 
[551 156j . An interesting application of a s+ii?s-series with the nontrivial power 
variable y appeared recently in [57] . 



5. An elliptic analogue of the Gauss hypergeometric function 

5.1. Definition of tiie ^-function and a connection v^ritii the root system 
El. 

The Euler integral representation for the 2^1(0, ^; c; a;)-function differs from the 
beta integral by the presence in the integrand of an additional term depending on 
two new parameters \\\ . An elliptic analogue of the Gauss hypergeometric function, 
which we shall be denoting by the symbol y(ti, . . . , ts;^, 9), is also given by a two- 
parameter extension of the elliptic beta integral |31| : 

y(ti,...,t8;p,9) = K / ^ , (5.1) 

where eight parameters ii, . . . , ig S C and two basic variables p,q € C satisfy the 
constraints \tj\, \p\, \q\ < 1 and the balancing condition 0^=1 tj — p^q^ ■ For other 
values of the parameters tj the ^-function is defined by analytical continuation of 
the integral (|5.ip . This continuation is build by the replacement of the contour of 
integration T by a contour C separating the sequences of poles of the integrand 
at z = tjp°-p^ , a, 5 = 0, 1, . . ., converging to zero from the poles at z = tj^p^'^p^'' , 
diverging to infinity; this does not assume now the restrictions \tj\ < 1. Shrinking 
the contour C, one can represent the resulting function as an integral over T and 
the sum of residues crossed by the contour during this deformation. 
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Suppose that a pair of parameters satisfies the condition tjtk = pq. Then the 
function is reduced to the elhptic beta integral (|3.ip . which follows from the 
reflection formula for the elliptic gamma function T(z;p,q). 

Let us consider symmetries of the 1^-function. Evidently, it is invariant with 
respect to the permutation of p and q and S's-group of permutations of tj related 
to the root system A-^. It appears that there exists a transformation extending S'g 
to the Weyl group for the exceptional root system Et. It is derived with the help 
of the double integral 

UU b,w^'-^ g) Tjcz^^w^^; p, q) dz dw 

T2 T{z^^,w^^;p,q) z w ' 

where Cj, bj,c £ C, \aj\, \bj\, \c\ < 1, and 11^=1 % ^ '^^ 11^=1 — PI- Computa- 
tion of the integrals over z or over w in different orders with the help of formula 
p.ip yields the fundamentally important relation [TTj 

V{t)^ n ^itjtk,tj+4tk+4;p,q)V{s), (5.2) 

l<j<fc<4 

where V{t) = V{ti, . . . ,t%;p,q) and 



Sj^etj, j = 1,2,3,4 ^ _ / pq _ t^tetjts 



Sj ^ £ tj, j = 5,6,7,8 ' y tit2t3t4 V pq 

and \tj\, \sj\ < 1. 

Repetition of transformation (|5.2p with the parameters 33^4.5^6 playing the role 
of ti,2,3,4 and subsequent permutation of the parameters t^jt^ with t^jt^ result in 
the equality 

4 

Vit)= llTit,tk+4;p,q)V{Ti/h,...,Ti/U,ui/t5,...,ui/ts), (5.3) 

where T = ^1^2*3^4, U = tstetrh and |T|i/2 < \tj\ < 1, \U\^/'^ < \tj+4\ < 1, j = 
1, 2, 3, 4. Equating the right-hand sides of equalities (|5.2p and (|5.3p and expressing 
parameters tj in terms of Sj, we obtain the third transformation 

V{s)= II r(s,5fe;p,g)-\/(v^/si,...,Vw/s8), (5.4) 

l<j<k<S 

where < \sj\ < 1 for all j. 

Let us connect parameters of the function V{t) with the coordinates of an Eu- 
clidean space Xj e by the relations tj = e^"^^^^ (pq)'^/^ . Denote as (x, y) the 
scalar product in and as - an orthonormal basis, {ei,ej) = Sij. The root 
system consists of the vectors v = {ei — Cj, i 7^ j}, and its Weyl group 5*8 - from 
the reflections x — > Sy{x) — x~2v{v,x)/{v,v), acting in the hyperplane orthogonal 
to the vector Y!h=i ^^^^ case the coordinates of the vectors x = 

satisfy the relation Y^^=i -^i ~ 0, which is guaranteed by the balancing condition. 

Transformation of the coordinates in (|5.2p corresponds to the reflection 5.„(a;) 
with respect to the vector v = (X]f=5 ~ '3i)/2 of the length {v, v) = 2. It 

extends the group Ss to the Weyl group for the exceptional root system Et. A 
different proof of equalities (|5.2p - (|5.4p is given by Rains in f34j, where a relation to 
the i?7-group is indicated for the first time. 
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In [TT] it was shown that the function is reduced to the product of two nVw- 
series for special discrete values of one of the parameters. Let us denote A = io • • ■ ^4; 
then 

y(to, . . . , is, ts" V+''z'"+\ V->'-'";p, 9) 

X 12 1^11 — ; — ,toti,tot2,tot3,tot4:,q ,— — ;q,p 
\ q q h , 

f Ato Au _„ Ap" \ , , 

X 12V11 ; ,toti,tot2,tot3,toti,p ,— — ;p, g , (5.5) 

\ P P h J 

where the contour of integration in the definition of ^-function is chosen in such 
a way that it separates sequences of the integrand poles converging to zero and 
diverging to infinity. For m — 0, we obtain an integral representation of a sep- 
arate terminating i2Vii-series. Since the left-hand side of (j5.5p is symmetric in 
to, . . . ,t4, the same should hold for the right-hand side as well, which follows from 
the symmetry transformation for series to be considered below. 

5.2. The elliptic hypergeometric equation. 

The addition formula for theta functions (|A.5p . being written in the form 

ts0{t7t^\tGZ^^;p) + te9{ts4\trz^^;p) + tre{tet^\tsz^^;p) = 0, 

leads to the connection formula 

teViqte) , tjV{qtr) , tsV{qts) 



0{tG4\te4^-p) e{trt^\t7tf^;p) e{tst^\tstf^;p) 



0, (5.6) 



where V{qtj) denotes the T^-function with the parameter tj replaced by qtj (this 
leads to the balancing condition 11^=1 *i — P^g)- Indeed, it is easy to check that 
the same equality holds for the ^-function kernel, after integration of which over z 
one obtains (|5.6p . 

A substitution of transformation (|5.4p in (|5.6p yields 

t&9{tT/tfi,tfi/U;p) tr9{te/t7,ts/tT,p) 
n-=ie(i8i,/g;p) 



tsOite/ts^tj/ts^p) 



Viq-Hs) = 0, (5.7) 



where 11^=1 *i = P^q^- 

Let us consider three equations appearing after the replacements tg qtg, in 
(|5.7p and q^^t^ or — > g^^ty in (|5.6p . Excluding from them the functions 

V{q~^tQ,qt^) and V{q~'^t-j,qt^), we obtain the elliptic hypergeometric equation 

A[h, ...,te, tr, ts, q;p) (u{qte,q-%) - C/(i)) (5.8) 
+A{h, ...,tr, te,ts,q;p) (u{q-he, qtj) - C/(t)) + U{t) = 0, 
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where 

A{ti,...,ts,q;p) . I . , . I vll~o777 T (^-^^ 

and 

V{t;p, q) 



Uit) :-{7a;p,g) 



uL6ntktr;p,q) 

The potential A{ti, . . . ,ts,q;p) is an eUiptic function of all its parameters, i.e. it 
does not change under the transformations tj — )■ ptj, j = 1,...,7, provided we 
count as a parameter depending on others through the balancing condition. 
Because of the symmetry of the [/-function in p and q, it satisfies another elliptic 
hypergeometric equation, which is obtained from (|5.8p by permutation of p and q. 
A change of variables brings the function A{ti, . . . ,tg,q;p) to the completely 5*8- 
symmetric form (j4.6p with arbitrary ui, . . . , ug and = 1 (a remark from D. Zagier 
to the author). However, the elliptic hypergeometric equation itself is maximally 
Sg-symmetric [5T1 [55] . 

Let us denote — cx, = c/x and pass to a new set of parameters 

efc = — , fc = 1, . . . , 5, £8 = —, £7 = —. 

ctk h q 

We fix the parameter eg from the requirement that the balancing condition takes 
the form nfe=i ^fc — P^Q^, which yields c = ^Je^eijp^. Now, after the replace- 
ment of C/(t) by some unknown function j{x) in ()5.8p . we can rewrite this elliptic 
hypergeometric equation in the form of a (/-difference equation of the second order: 

A{x) ifiqx) - fix)) + A(x-i) {f{q-'x) - f{x)) + vf{x) = 0, (5.10) 

Ai^) = ^F^^, -n^(— (5-11) 

e{x\qx^]p) \ q J 

where one can explicitly see the Sg-group of symmetries in parameters. We have 
already one functional solution of this equation: 

V{q/ c£i, . . . ,q/cs^,cx,c/x,c/£g ,;p, q) 
T{c^x^'^ /e%,x^'^ei]p, q) 

In order to build other linearly independent solutions, one can use symmetries 
of the equation (|5.10p which do not represent symmetries of the function (|5.12p . 
For instance, the second solution can be obtained by multiplication of one of the 
parameters ei, . . . , £5 or x by the powers of p or by permutations of £1, . . . , £5 with 
£g. 

Since the function A{t^ q;p) in equation (|5.8p does not change after the replace- 
ments tg p^^tQ^t-; pti, the function U{p^^tQ,pt7) also represents a solution 
of this equation. The Casoratian (a discrete Wronskian) of these two solutions was 
computed in the paper [SS]. Let us multiply equation (|5.8p by U{p~^tQ,pt-j), and 
the equation for U{p~^te,ptY) by U{t) and subtract one of the resulting relations 
from another. As a result, we obtain the equality 

A{ti, . . . ,te,tT,ts,q;p)D{p^^te,q^^tr) = A{ti, . . . , ty, ig, tg, 9;p)D(p"^q~4g, ir), 

(5.13) 

where D denotes the Casoratian 

-0(^6, t?) = U(pqt6,t7)U{te,pqt7) - U{qte,pt7)U{pt6,qt7). 



hix;e;p,q) = t...2^±i ^±1. -^ ' (5-12) 
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After substitutions tg — > ptg, tj — > qtj the balancing condition takes the form 
tail — pq/ti ■ ■ -tsts, and ()5.13p can be considered as a q-difference equation of 
the first order in ty. It is easily solved, and the solution is defined up to the 
multiplication by an arbitrary g-periodic function ip{qt-/) = (/'(^t)- 

Repeating the same procedure with the equations appearing after the permu- 
tation of the parameters p and q, we find (^(pty) = ^{tr), so that 1^9(^7) does not 
depend on t^. Further investigation of the structure of pole residues of the V- 
functions in D{t(,,t'i)^ which are crossed by the integration contour in the limit 
t-j tg^, fixes completely the form of ip and leads to the identity [59] 

Vipqte,tr)V{U,pqt7)-t,%^V{qU,pt7)V{pte,qt7) = '^^ 

(5.14) 

where and ^7 can be replaced by any other pair of parameters. 

The described solutions of the elliptic hypergeometric equation exist for \q\ < 1, 
whereas the equation itself (|5.8p does not demand such a restriction. Because of 
the symmetry 

/pi/2 pl/2 \ 

A -j^,q;p] ^A{ti,...,ts,q ■,p) , 

the transformation tj p^^-^/tj, j — 1,...,8, leads to the change of the base 
variable q ^ q~^ in (|5.8p . This yields the following solution of the elliptic hyper- 
geometric equation for |g| > 1 [59| 

Ilk=6^ip/'^kts,ts/tk;p,q ^) 

In order to build solutions of equation (|5.8p (or (|5.10p ) for \q\ — 1, one can 
use the parameterization of the base variables (|2.6p . (|2.7p . make substitutions 
tk — e^^'f*/'^^ and repeat the whole chain of arguments given above with the 
replacement of T{z;p,q) by the modified elliptic gamma function G{u;uj). In the 
same way as in the case of modified elliptic beta integral, we obtain again the 
^-function, but with a different parameterization. This procedure appears to be 
equivalent to the application of the modular transformation (0^2, ^^^3) ~* (—^3,^2) 
to solutions described above. Indeed, the potential ^(e2''*fl/"^ . . . , e^''*S8/'^2 ^ g.^) 
is invariant under this transformation, but U{ti, . . . ,ts]p,q) gets transformed to 
jj^g-2Trigi/iJ3^ . . . , e^^TTigg/ws.^^ fiy 'pjjg Jitter function provides thus a new solution 
of the elliptic hypergeometric equation well defined for jql = 1. It is evident that 
this function satisfies also a partner of equation (|5.8p obtained from it by the per- 
mutation of ui and uj2 ■ An example of a similar situation with two equations for 
one function at the g-hypergeometric level is given in [60] . 

Different formal degenerations of the 1^-function to g-hypergeometric integrals 
of the Mellin-Barnes or Euler type are briefly considered in [31] [55] . A detailed 
and rigorous analysis of the degeneration procedure is performed in |281 161j . It is 
necessary to note ^58) that the elliptic hypergeometric equation emerges as a partic- 
ular case of the eigenvalue problem equation for the one particle Hamiltonian of the 
quantum model proposed by van Diejen [62| and investigated in detail by Komori 
and Hikami [63] . This Calogero-Sutherland type model represents a generalization 
of the Ruijsenaars [53] and Inozemtsev [5S] systems. 
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6. Chains of symmetry transformations for functions 

Symmetry transformations for the plain and g-hypergeometric series are built 
with the help of the Bailey chains [TJ [SB] . This technique was generalized to the 
elliptic level in [52l [67] . Let us sketch this generalization. 

Two sequences of numbers Q!„(a, /c) and /3„(a, /c) by definition form an elliptic 
Bailey pair with respect to the parameters a and k, if 

f]n{a,k)^ ^ Mnm{a,k)am{a,k), (6.1) 

0<m<n 

where 

j,j _ 6>(fc/a)„-,„6l(fc)„+„, 6l(ag^";p) 

In the matrix form [3{a,k) — M{a,k)a{a,k), where a and /3 denote the columns 
formed by a„ and /?„. 

Let us introduce the diagonal matrix 

Dnm{a; b, c) = Dmia; 6, c)(5„„i, An(a; c) = ^^^^j^' ^^J^^,^ (^^) . (6.3) 

Theorem 4. Lef a{a, t) and P{a, t) form an elliptic Bailey pair with respect to the 
parameters a and t. Then the quantities 

a'{a,k) = D{a;b,c)a{a,t), f3'{a,k) = K{t,k,b,c)P{a,t), (6.4) 
K{t, k, b, c) := D{k; qt/b, qt/c)M{t, k)D{t\ b, c), 

where qat = kbc and b, c are two arbitrary new parameters, form a new elliptic 
Bailey pair with respect to a and k. 



Proof. Substitution of (j6.4[) and of the relation /? = Ma in the required equality 
P' = Ma' leads to the matrix identity 

M{a, k)D{a; b, c)M {t, a) = D{k; qt/b, qt/c)M {t, k)D{t; b, c). (6.5) 

After substitution of the explicit expressions for matrices, one can see that it is 
equivalent to the Frenkel-Turaev summation formula (|3.8p . □ 



Because M„m(a, a) ~ Snm and Dnm{bc/q; b, c) = (5„,n, we find after setting t — k 
in ()6.5|) that M(a, k)M{k, a) = 1. The inversion of the matrix M is reached thus by 
the permutation of parameters a and k (in the p = case this fact was established 
in |68| : a more detailed discussion of such matrix inversions is given in [53, 56 , 69J). 
Therefore, a{a, k) — f3{k, a) and /3(a, fc) = a{k, a) define new Bailey pairs to which 
one can apply the transformation (|6.4[) . The described rules of composition of new 
Bailey pairs generate a binary tree of identities for different products of matrices 
M and D, which are equivalent to some nontrivial identities for (multiple) elliptic 
hypergeometric series. 

From the relation M{a, f)M{ f, a) — I we find the simplest Bailey pairs an '' (a, /) = 
Mniif, a) and (3n^ (a, /) = 5ni- Let us set a' {a, t) = D{a; d, e)M{f, a) and /3'(a, t) — 
K{f,t,d,e), where qaf = tde and it is assumed that a' and p' are the matrices 
whose columns form the Bailey pairs. Then the equality P'{a,t) = M{a,t)a' {a,t) 
is equivalent to (|6.5p . The relation f3"{a,k) — M{a,k)a" {a,k) with a"{a,k) — 
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D{a;b,c)a'{a,t) and f3"{a,k) — K{t,k,b,c)P'{a,t), where qat = kbc, leads to the 
identity 

^^ t h ^ I. n -n f-i N 0{qa,k/t,qt/b,qt/c)n 
i2Vii{a; b, c, d, e, kq ,q ,af ■,q,p)- 



9(k/a, qt, kb/t, kc/t)n 
X uVnit; b, c, td/a, te/a, kq"", q-'\tf-^;q,p), (6.6) 

where / — kbcde / (aq)'^ . This relation represents an elliptic analogue of the Bailey 
transformation for terminating io</39-series which was proved for the first time 
in [S] by a different method. There exists also a four term Bailey transformation for 
non-terminating loV'g-series [3]. Its elliptic generalization is given by the ^-function 
transformation ()5.2|) (written in the integral form since its infinite series version is 
not well defined). 

Integral analogues of the Bailey chains were discovered in the paper [30| . where 
a number of symmetry transformations for the elliptic hypergeometric integrals has 
been built with the help of this technique. The functions a{z,t) and P{z,t) form 
by definition an integral elliptic Bailey pair with respect to the parameter t, if they 
are connected to each other by the relation 

(3{w,t)^K / T{tw^'^z^'^;p,q)a{z,t) — . (6.7) 
Jf z 

Theorem 5. Let a{z, t) and P{z, t) form an integral elliptic Bailey pair with respect 
to the parameter t, |i| < 1. We take the parameters w, u, s satisfying the conditions 
w €T and \s\, \u\ < 1, \pq\ < \t'^s'^u\. Then the functions 



, T{tuw^^;p,q) 



, T{t^s'^,t^suw^'^;p,q) [ T{sw^^x^^ ,ux^'^-p,q) dx 
l3{w,st) = K ^^^^^ „ ^ , / „^ Pi^^^y 



T{s'^ ,t'^ , suw^^\p,q) Jf T{x^'^,t'^s'^ux^^;p,q) x 

define a new integral Bailey pair with respect to the parameter st, and the functions 

,, , Tis^t^ ^uw'^^\p,q) f Tit^sux^^ ,sw'^^x^^p,q) , .dx 

a {w,t) ^ K — — / — — a[x,st} — , 

L (s'',t^,w^'',t''s''uw^'-;p,q) Jf 1 (sux^'^;p,q) x 

^("'^)^ r(t.Wi;^,,) ^^"'^^^ 

define a new integral Bailey pair with respect to the parameter t. 

In order to prove the first statement, it is sufficient to substitute relation (|6.7p 
in the definition of (3'{w,st), to change the order of integrations, and to use the 
elliptic beta integral (|3.ip . The second statement is proved in an analogous way. 

In the paper [3S] it is shown that under certain restrictions the integral transfor- 
mation a{x,t) — + P{x,t) has a very simple inversion. Let p,q,t € C are such that 
\p\, \q\ < |tp < 1. For a fixed w £ T, we denote as a contour inside the annulus 
A = {z E €.\ \t\ ~ e < \z\ < \t\^^ + e} for some infinitesimally small positive e, such 
that the points t^^w^^ are lying inside Cw Let f{z,t) be a holomorphic function 
in A satisfying f{z,t) = f{z~^,t). We define an integral transformation 

f dz 
g{w,t)^K S{z,w;t-^)f{z,t) — , (6.8) 
.In. Z 
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where 

Then for \t\ < \z\ < we have 

f d 
f{z,t) = K / S{w,z;t)g{w,t) — . (6.9) 

This relation coincides in essence with the definition of eUiptic integral Bailey pairs. 
It can be shown that the two ways of building chains of integral Bailey pairs indi- 
cated above are related to each other by the described inversion of integral trans- 
formation ()6.8|) . 

With the help of Theorem [5] one can build infinite sequences of Bailey pairs 
starting from a given initial pair. The simplest pair can be built with the help 
of the elliptic beta integral p.ip . Each new application of the substitutions indi- 
cated above brings in two new parameters. Equality (|6.7p . being applied to the 
appearing new Bailey pairs, leads to a binary tree of identities for multiple elliptic 
hypergeometric integrals with many parameters. 

As an illustration, we describe a chain of nontrivial relations for the integrals 

T(rn)u , -s f Ilj=l^'^{tjZ^^;p,q) dz ^J^^, I xm+i 

1^ ■ ■ • ,<2m+6) = K / ^, ^2 T , M = W ^ , 6.10) 

Jt r(z±^;p,g) z jj-^ 

where \tj\ < 1. With the help of the elliptic beta integral (|3.ip it is easy to verify 
validity of the recursion 

. n 2m+5<k<l<2m+S 

T{tkti;p,q) 

r(p„i;p, q) 

XK/ ^, , „ r 1"- '[ti,...,t2m.+4:,PmW,PmW ) , 

where = Ylk^tm+b^k/pq- This equality gives a concrete realization of the 
Bailey pairs a ^ /'™^ and (3 ~ after a change of parameters. For m = 

0, substitution of the explicit expression for l'^*'^ (|3.ip to the right-hand side of 
(|6.1ip leads to the identity (|5.2p . Another important consequence of recursion 
(|6.1ip is considered in the next section. In general equality (|6.1ip yields an m- 
tuple integral representation for J^™) analogous to the Euler representation for the 
m+iJ^m-function. 

7. BlORTHOGONAL FUNCTIONS OF THE HYPERGEOMETRIC TYPE 

7.1. Discrete biorthogonal functions with the continuous measure. 

Let us denote £{t) :— i2Vii(tQ; io<i, . . . ,totj;q,p) with the balancing condition 
nL=o ^™ ~ 1^ ^^'-^ ^^^^ termination condition to^fc = 9"" for some k. The contigu- 
ous with £{t) functions, which are obtained by the change of parameters ti and tk 
to tiq and t^q^^, will be denoted as £{tf,t^). Then, using the addition formula 
for theta functions, it is not difficult to check the equality [21 SS] 

_ S{qtl,qHl,qt7/tQ,tQtT/q]p) A _^Mrjp) ^^ 

11 flr«^„/f„-r,^^^*0.t6j- U-iJ 



0{qto/t6,qHo/te,to/t7,t7/qto;p) ^J-^ 9{qto/tr;p) 
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Substitution of an elliptic analogue of the Bailey transformation and of its iterations 
in ()7.1|) allows one to build many formulae of such a type. One of them has the 
form 



e{tG/qto,te/qHo,te/t7;p) ^J-^ 



= %.§,^^.g;p) n%^"A-;^)^(t)- (7-2) 

Let us replace in (|7.ip the parameter Iq by gig, and ti by tj/q, and substitute 
£{tQ ,t^) from the resulting equahty and £{1^,1^) from (|7.ip in (|7.2p . This leads 
to the equation 



e{qtr/te,tr/te;p) 



+ 0(gA6i7;p)n^(^otr;p)f(t) = 0. (7.3) 



5 



These relations are analogues of equations (|5.6p , (|5.7p , and (|5.8p for the ^-function, 
and they can be obtained from them by the residue analysis for some sequences of 
poles of the integrand. Let us replace in (|7.3p the ^-function by 

D / N i ^ 9 9 -„ Aq""-^ £6 \ .X 

i?„(a;;g,p) = i2Kii — ; , , ,g , ,—,£(,x;q,p\, (7.4) 

yes EiEs e2£8 e3£8 £8 / 

where n = 0, 1, . . . and A = e\e2Sze^e^, and denote q^" ~ pqle^e^, Aq"-^^/es — 
pq/e^Es- Then, after the change of notation for parameters = £6/£8, ■ • ■ , toh — 
EqX, we obtain equation (|5.10p with = Sg/q and discrete values of one of the pa- 
rameters. The function i?„(x;g,p) defines thus a particular solution of the elliptic 
hypergeometric equation obeying the property Rn{px; q,p) = Rn{x; q,p). (Nota- 
tions of the paper where this function was investigated, pass to ours after the 
changes io,i,2 ^ £1,2,3, ^3 £6, U - > £8, ^ SiSs/pq and A^i/qt^ pq/e^es.) 

The parameters ei, . . . ,£6 enter (|5.10p symmetrically. Because of the balanc- 
ing condition 0^=1^*; ~ P^Q^, the function Rn{x;q,p) is invariant with respect 
to changes Sk — > p£k, k — 1,...,5. This guarantees the symmetry of (|7.4p in 
£1, . . . , £5, and any of these parameters can be used for the series termination. Per- 
mutation of one of the parameters £1,2,3.5 with £g and application of the elliptic 
Bailey transformation for terminating i2VLi-series leads to Rn{x;q,p) up to some 
multiplier which does not depend on x. 

We identify now the parameters in equation (|7.3p in a different way: 

to — — , ^0*1,2,3 — — , tot^ — — , io^s — ^6^1 hte = 9 totj ~ — 

£8 £i,2,3£8 a; £8 
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This leads to a three term recurrence relation in the index n. 



n^(er£6g";p) \p] {Rn+l-Rn) 



esq 



-;p 



\ /^2-n^itl 



£8 ' £8 ^ ^ j-r ErSeesq^ \ a 1 \fu o\ 



{'-^,e,x^';p) n (^■,p)Rn = (7.5) 



with the initial conditions = and i?o = 1. 
Let us introduce the functions 

where ^ and are arbitrary gauge parameters, f 7^ -q^^p^ , A; G Z. Then relation 
(|7.5p can be rewritten in a more structured form 

(z(a;) - a„+i)p(yl(7""Ve8) (^n+i(a;; - -R„(a;; 

+ (z(x) - /3„_i)p(g"") (i?„_i(a;;g,p) - i?„(x;q,p)) 

+ <5(z(x)-z(£6))i?n(a:;g,p) = 0, (7.6) 



where 



^ 9£6. qt qt qt 

' £8* ' £8* ' £l£2 ' £2£3 ' £l£3 ' 



9*^ £8 2i*££8 . „ 



q^£6 q 



,£6V'^^;P 



A EiEs £2^8 £3£8 

The initial conditions — and Rq — I guarantee that Rn{x;q,p) are rational 
functions of z{x) with the poles at the points ai, . . . ,«„ (i.e., all the dependence 
on a; enters i?„ only through the variable 2(2^))- 

Suppose that 4'\ is a solution of an abstract generalized eigenvalue problem 
Vicpx = X'D2(j)\ for some operators I?i,2- Let also a scalar product is given, 

which defines the formal conjugated operators 2?^2 t>y the standard rule {T>l2'^\<t>) = 
(i/)|2?i_2(^). Let ip\ denote solutions of the dual generahzed eigenvalue problem 
P^Va' = APJVa. Then = (V'^lPi " A2?2)0a) = (m - A)(I?i^i/;^|0A), that is 
the function 2?JV'm orthogonal to (/>a for /i 7^ A. Consequences of this well 
known fact of the linear algebra were investigated in detail by Zhedanov [70j in 
the case when I?i^2 are the Jacobi matrices (i.e., the tridiagonal matrices). In par- 
ticular, such generalized eigenvalue problems were shown to be equivalent to the 
theory of biorthogonal rational functions generalizing the orthogonal polynomials. 
They are connected also to the recurrence relation of the i?77-type investigated 
in [71j and to the orthogonality relations appearing within the theory of multi- 
point Pade approximation [72l |73] . Recurrence relation (|7.6p belongs to this class 
of problems and, therefore, there exists a linear functional C with the condition 
i^{Tm{x;q,p)Rn{x]q,p)} — hnSnm for some rational functions Tm and normaliza- 
tion constants hn- 
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The elliptic hypergeometric equation for i?„-functions can be rewritten in the 
form of a generalized eigenvalues problem 

'D{e'^,e'^)Rn = \n'D{e'l,e'l)Rn, 

where 

I?(£4,e5) - A{x){T,,^ - 1) + A{x-^){T-^l - 1) + i^, T,,,f{x) = f{qx), 

denotes the operator permitting to rewrite equation (|5.10|) as I?(e4, e5)/(a;) = 0. 
The primed parameters are arbitrary under the restriction £465 — e\e'l = £465 , and 
other parameters remain untouched (that is why the dependence on them is not 
indicated). The spectral variable 

^ 6>(£4/g4,£4/£5;P) 
" 0(£4M',£4/4';P) 

is discrete because £4 = •pq^^'^ je^. If we take as the functional £ the integral whose 
kernel coincides with the kernel of elliptic beta integral, 

£{0(2;)V^(2;)} = K /" n^-i 2,3,6,8 (^J 

JT i (,2^ J ,p, 9J X 

then the functions 

rr ( \ T/ (^^^ AAA £6 \ 

r„(x;g,p) = i2Vii ; — , — , — ,£52;, — ,9 , ;<7,P , (7.7) 

\ g £1 £2 £3 X £8 / 

serve as an analogue of PJV'm fo'^ (-^5 9' P) • 9' ^^^^ the rational functions of 

z{x) with the poles at the points (3i, . . . , Pn which are obtained from Rn{x\ q,p) after 
the parameter change £§ — ^ Pi/ A (dependence on p in the parameters disappears 
because of the ellipticity of the 12 Pi 1 -series in them). 

We denote i?„m(x) := Rn{x;q,p)R^{x;p,q) andT„„(a;) := Tn{x;q,p)T^{x;p,q), 
where all 12 Vii -series terminate simultaneously because of the modified termina- 
tion condition £4£8 = ^"1+1^,"+!^ n, m = 0, 1, . . .. Since Rmiqx;p,q) — Rm{x]p,q), 
the functions Rnm represent now solutions of not one, but two generalized eigen- 
value problems differing from each other by permutation of p and q. Therefore 
the orthogonality relations in our case appear to be more complicated than for the 
biorthogonal rational functions. 

Theorem 6. The Junctions Rnm{x) and Tnm{x) satisfy the following two-index 
biorthogonality relations: 

Tnl {x)Rrnk (x) ^^±2 ^ j^^±l . " = ^^1 , (7.8) 

where S ~ {1, 2, 3, 6, 8}, Cmn,ki denotes a contour separating the sequences of points 
ejp''q''{j = 1,2,3,6), £8p"~'=g''~",p"+i"'g''+i~"M, a, b ^ 0,1, from their x ^ 
x~^ reciprocals, and the normalization constants have the form 

Ki = nj<fc,j,fcgs '"'^'^^ h^{q,p) ■hi{p,q), 
Y[j^s^{Ae^ ■p,q) 

, , . e{A/qes;p)9{q,qe(i/es,eie2,sie3,e2e3,AeG)nq~'' 
hn{q,p) - 



9{Aq'^'^/q£s;p)9{l/seSs, £i£6, £2£6, ^s^e, A/qse, A/qes)n ' 
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A direct proof of this statement by a straightforward computation of the integral 
on the left-hand side with the help of formula (13. ip and the Frenkel-Turaev sum 
is given in [11]. Appearance of the two- index orthogonality relations for univariate 
functions is a new phenomenon in the theory of special functions. It is worth of 
noting that Rnm (x) and T„m (x) are meromorphic functions of x e C* with essential 
singularities at x = 0, oo; only for fc=:/ = Oorn = m = they become rational 
functions of some argument depending on x. For k — I ^ 0, one can take the limit 
p — > with fixed parameters and obtain the functions Rn{x; q,0) and T„(x;g, 0), 
which coincide with the family of continuous iQifg biorthogonal rational functions of 
Rahman [39^ . A further degeneration of these functions leads to the Askey- Wilson 
polynomials [41j . Additional restrictions on one of the parameters in Rn{x; q,p) and 
Tn{x; q,p) leads to a finite-dimensional systems of biorthogonal rational functions of 
a discrete argument [9l [T^ , generalizing the Wilson functions ^75^ . An elementary 
approach to the analysis of these functions, related to the elliptic 6j-symbols [5], 
was suggested by Rosengren in |76| . Some properties of the functions Rn{x; q,p) 
are investigated in the recent paper [48j . 

One can build a relation analogous to (|7.8p on the basis of the modified el- 
liptic beta integral (|3.10p [31] . For this it is necessary to use the parameteriza- 
tion of base variables in terms of the quasiperiods ^^-"1,2,3 and pass to the func- 
tions r„(u; wi,u;2, W3) = i?„(e2""M . g277ia)i/i^2 ^ £277*0)3^)^ where we have also sub- 
stituted Ej — ^'"^ail^'i^ Analogously, it is necessary to redenote Tn{x;q,p) as 
Sn{u;uji,uj2,uj3) a.ndhn{q,p) as /i„(wi, a;2, ^3)- The products r„„ (u) = r„(u; wi, 0^2, W3) 
rm(u;a;2, wi, W3) and s„m(u) = s„(m; wi, 0^2, a;3)sm(u; 072, t^i, ^3) are invariant with 
respect to the permutations lui lu2, n ^ m. Then, for a specially chosen contour 
ki, we have 

n,g5 G{gj ± u; uj) du 

Snl{u)rrak{u)-4--— ^— r = hnl Smn hi, (7.9) 

c„„.,fc, G{±2u,A±u\u)) U2 

where A — J2j(ES 9j 

7 nj<m,j.mes^(ffi +5™;^) , , , . 

In distinction from the previous case, the limiting transition to the g-hypergeometric 
level Im(cj3/a;2), Ini(cj3/a;i) -~^ +00 (that is, p,r —^ 0) is well defined and preserves 
the two-index structure of biorthogonality relations. In particular, the rnm{u)- 
function degenerates now to the product of two q-hypergeometric series 

r„„(u;u;i,c^2) = ioW'9(e2"*(^«-3*^)/'^^e2"('^^-^^-»«)/"% . . . ,e2-(^«+")/^ 

X -^Qj4/g^g2'ri(96-ffs)/i^l.g27ri(w2-3i-g8)/wi^ _ _ ^ ^ g27ri(c/6+u)/(^i . ~-l ^ 

whose basic variables are related by a modular transformation and the normaliza- 
tion of the measure is given by the integral (|3.13p . 

7.2. A terminating continued fraction. 

A terminating continued fraction related to the rational functions Rn{x; q,p) 
is computed in the paper [IB]. Let J7„ and Vn denote two sequences of numbers 
satisfying the three term recurrence relation 

V'n-l-l = CnV'n + ?7nV'n-l, n=l,2, ... (7.10) 
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with some coefBcients ^„ and and the initial conditions Uq ^ 0, Ui = 1 and 
Vq = 1, Vi = ^0- It is well known that their ratio is related to the finite continued 
fraction 

^ = ^ , n=l,2,.... (7.11) 

Vn , 7?1 

4o 



£.1+ ■■■ Vn-l 

Let us define polynomials of z{x) of the n-th degree: 

n n — 1 

F„(z(x)) = K„ - ak) Rn{x;q,p), k„ = J| p{aq^^^) 

k=l j=0 

and set Po{z{x)) = 1. Replacing Rn{x;q,p) by P„(z(a;)) in (|7.6p . we obtain the 
following recurrence relation 

F„+i(z) + {v„ - Pnz)Pn{z) + Un{z - a„)(z - /3„_ i )F„_ i (z) = (7.12) 

with the initial conditions P_i = 0, Pq = 1 and the recurrence coefficients 

Un - p(g-")p(^g"-Ve8), Pn - p(^g"-Ve8) + ^(9"") " ^, 

t'n = a„+i/9(Ag"-V£8) + I3n^ip{q-'') - 5z{eQ). (7.13) 

In this case Vn — Pniz), and J7„ = P^^^j^(z) are the associated polynomials of the 
degree n — 1 in z. 

Let us suppose that the polynomial Pn+i{z) has only simple zeros, that is 
Pn+i{z,) = 0, Zs = zi^+^\ s = 0,1,..., TV, z, ^ Zs' for s ^ s' . Then the 
corresponding continued fraction can be expanded into the partial fraction (as a 
rational function of z) 

£nA^ = V qs = ^ ^ (7 14) 

The Casoratian of any two solutions C/„ and Vn of (|7.10p satisfies the relation 

Un+lVn - UnVn+1 = (-l)"r,i • • • r^n{UlVo ~ UoVi) , 

which yields 

P„(z)P^l)(z) - P„+i(z)PiL\(z) = hnAn{z)Bniz), (7.15) 

where hn — uiU2 • • - u^ and 

n n 

An{z)^\{{z-a,), P„(z) = [](z-ft_i). 

■i=l i=l 

Fixing n = N and setting z = z^, s = 0, 1, . . . , iV, in (|7.15p . we can express pIJ\zs) 
in terms of PNizs), h^, An{zs) and Bm{zs)- This leads to the following convenient 
for computations expression for the residues of the poles gs in (|7.14p : 

hNAN{Zs)BN{Zs) 
P'j^j^^{Zs)Pn{Zs) 

In the limit 

e3e6 = 'Z"'^+% iV = 0,l,..., e^O, (7.17) 
we find that un+i — >■ 0, and the continued fraction terminates automatically. It 
appears that this fraction can be computed in the closed form using formula (|7.14p . 
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For £ 0, the rational function Ri\i+i{x; q,p) diverges, since the elhptic Pochham- 
mer symbol 0{s3e6)N+i — » in the denominator of the last term of the i2Vii-series. 
However, we have simultaneously kat+i 0, so that the polynomial Pn+i{z) takes 
the finite value, and its zeros are found explicitly: Zs = z{seq'^), s = 0,.. . ,N. It 
appears also that the polynomial Pn{zs) is computable in the closed form owing to 
the Frenkel-Turaev summation formula. The other quantities defining the residues 
gs are found sufficiently easily, although they are given by rather cumbersome ex- 
pressions. 

Suppose that the conditions of the simplicity of zeros Zs are satisfied as well 
as other restrictions on parameters guaranteeing that Uk ^ Q for k — I, . . . , N, 
the descriptions of which we skip. Then the terminating elliptic hypergeometric 
continued fraction has the following explicit representation: 



ui{z - ai){z - (3o) 

poz - vq 

Piz-vi-... UN{z-aN){z- Pn-i) 

PNZ - VN 

1 ( gee qx q _^ e6q^+'^ \ 

~ T7 — ^ r w: 12^^11 — q, se^i, £e£2, — , — ,q , ;q,p , 

where in the expressions for all recurrence coefficients, including ^, it is necessary 
to substitute £3 = q~^ js^ and z{x) — 9{x£,'^-^;p)/9{xri'^^;p). In 46 this result was 
presented in the different (additive) system of notation. 

This formula describes the most general terminating continued fraction of the 
hypergeometric type, which was found to the present moment. For fixed values of 
parameters, in the limit p one obtains the terminating continued fraction of 
Gupta and Masson [77] (see Corollary 3.3) described by a very-well poised balanced 
iQi/'g-series. Further specification of parameters leads to the continued fraction of 
Watson which, in its turn, is a g-analogue of the famous Ramanujan continued 
fraction (see the details in [78]). 



7.3. Continuous biorthogonality of the F-function. 

The relations described above (|7.8|) correspond to discrete values of one of the 
parameters in the elliptic hypergeometric equation. In paper [79], it was shown 
that the ^-function with general set of continuous parameters obeys also some 
biorthogonality relations characteristic to the continuous spectra. 

Let us consider the m = 2 case in recursion (|6.1ip . After imposing the con- 
straints on parameters t^tj = tgtg — pq the integral on the left-hand side becomes 
explicitly computable. Then, after a number of notational changes and application 
of transformation (|5.2p . there appears the equality 

f duj 
(/)(a;; c, dl^; s) = K / i?(c, d, a, 6; x, w|s)0(u'; a, 5|f ; s) — , (7.19) 

where the basis vectors have the form 

<j^{w, a, b\e, s) = T{sat\ sb^^\ &±i^±i;p, q) (7.2O) 

V ao 
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and 



i?(c, d, a, 6; x, w\s) 



ab 



cd 



cd 



as bs y pq 



pq 



Here the parameters a, b, c, d, s and x, ^ are arbitrary, but their choice should match 
with the condition that taken contours of integration separate converging to zero 
and diverging to infinity sequences of poles of the integrands. The variable ^ enters 
only the basis vectors 4'{w\ o, 6|^; s) and the kernel R can be considered as a "rota- 
tion matrix" with continuous indices x and w, which permits to change arbitrarily 
the parameters a and b. Equality (j7.2ip represents an integral generalization of 
the relation which was used by Rosengren in [76j for derivation of properties of the 
elliptic 6j-symbols. 

Denoting w = e*^, y = cos0, and using the equality 



Jt ™ Jo 



f{cos9)d9^2 J Jiy) 



dy 



we can write 



x; c, d|^; s) = 2iK / R{c,d,a,b;x,e^ \s)(f){e^ ;a, 6|^;s) 



dy 



y 



(7.22) 



In the limit c —f a and d b, there appears the following relation in the distribu- 
tional sense 



lim fl(c,d,a,6;e-^e"^)^ \ S{v-y), (7.23) 

where v = cos ip. 

A double application of relation (|7.19p with different parameters leads in an 
evident way to the self-reproducing property for the kernel 



K I R{a,b,c,d;x,'w\s)R{c,d,e,f;w,z\s) — = R(a,b, e, f; x, z\s) 

w 



(7.24) 



and the biorthogonality relation 



dy 



R{a,b,c,d]e"^,e"'\s)R{c,d,a,b;e"',e''^ \s)- 

1 vi-y^ 



(7.25) 

where v = cosi^j and v' = cosip'. Substitution of the expression for _R- function 
(|7.2ip in (|7.25[) results in the equality 



1 



pq 



pq 



pq pq / cd 



/ cd 



ir(e±2^e;p,g) 



V sa,sb,jq-e'^,jq-e~''^,i^,q-,J—e''',J—e 



ab 



ab 



c d 



X V \ sc. sd, 



[pq te [pq -te pq pq 



p 

I ab 



-lip 



dy 



cd 



as' bs' y pq 



[ab pq cd pq ±2^.p 
^pq ab pq cd 



(7.26) 
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The parameters v and v' can be considered as continuous spectral variables in the 
operator formulation of the elliptic hypergeometric equation. Therefore relation 
(|7.26p should follow from the latter equation, but the precise connection between 
them is not established yet. 



8. Connection with the Sklyanin algebra 

In [34j , Rains introduced an interesting finite difference operator connected with 
the root system BCn- For n = 1, it can be represented in the form 

0{az,bz,cz,dz;p) , 9{az-^,bz-^,cz-^,dz-^;p) -^^/^ 

(8.1) 

where T^q^'^f{z) — /(g^^/^z) is the g-shift operator and a,b,c,d are arbitrary 
parameters. Later on Rains noticed also [35l [80] that this operator is equivalent to 
the general linear combination of four generators of the Sklyanin algebra Sq, . . . , S3 

Defining relations of the Sklyanin algebra have the form 

SaSi3 — SpSa — i{SoSj + SjSq), 

SoSa — SaSo = iJpj{SpSj + SjSp), (8.2) 

where J^j are the structure constants of the algebra and (a,/?, 7) is an arbitrary 
cyclic permutation of the triple (1, 2,3). A representation of Sa as finite difference 
operators has been found in [82^ : 

Sa = *'°'^^7^ {9a+i{2u-2g\T)e^^^' ~ 0,+i(-2^ - 2.g|T)e-"^") , 
0i(2u\t) 

where e^^^" f{u) — f{u ± 77), and under the quantization condition g = irj, £ = 
0, 1/2, 1, . . . , their action in the space of theta functions of the order A£ was de- 
scribed. The combination of the generators 

n-=i^iK+a4 + 2.g)^ n-=i^iK+«4 + 25+i)^ 
2A[ai,a2,a3,a4) := ^-r-. "^o 7-^ — -^^^ 



3 



■/X n n ^ n,--1 ^1(0? + 04 + 2g + ^) 

with the normalization X]j=i '^j — "^y, can be represented in the form [80j 

A/ ^ nl=i^i(aj+w) a„ , 111=1 ^'i(aj-w) 

A ai,a2,Q3,a4 = a ^ + q f n \ ^ ' 

After the transition to multiplicative system of notation 

(a, 6, c, d) := e^''"!'^-^'-*, p ~ abed = e^^^^ z := e^^™, q := e""" 
there appears the operator described above (|8.ip : 

A(ai, a2, as, 04) = («P^^*(p;p)oo) e^''*^i:'(a, b, c, d;p; g). 
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The standard eigenvalue problem Dip — Xip appears to be very complicated, 
since it represents a difference analogue of the Heun equation |79| . On the one 
hand, it is known that the eigenvalue problem for the one particle Hamiltonian of 
the Inozemtsev model ^65J is equivalent to the Heun equation. On the other hand, 
the classical equations of motion for this Hamiltonian with the modular parameter 
T considered as a time variable leads to the Painleve VI equation [83]. The D- 
operator itself appears to be related to the van Diejen model [MllIS]. All this and 
the connection with elliptic hypergeometric functions described below demonstrate 
some mathematical universality of the operator (jS.ip . 

Consider the generalized eigenvalue problem of the form 

D{a, b, c, d;p; q)f{z; w; q^^^a, q^^^b; p) = A(w)D(a, 6, c', d'\p; q)f{z; w; q^^'^a, q'^^^b; p), 

. . (^-^^ 

where cd = c'd' . This equation is solved explicitly. Using the parameterization 

6{w^c/d, WyJ d/c\ p) 



Xiw) = , , 

Q(w^f7Jd! ,w^fdrid\py 

for |g| < 1 we obtain [79] 



/(z; a, 6; p) = r f ^z±\ ^-z^\ &±iz±i;p, ^ , (8.4) 



D {a,b,c,d;p;q) = -rj^D , — — , , —^;p]q 



up to the multiplication by an arbitrary function ip{z), >f{qz) — </'(z). Evidently, 
function (|8.4p coincides with the basis vector (|7.20p after a change of parameters. 
For z — > z^^ invariant functions, ipi^) = we define the scalar product 

Then the formally conjugated to D operator has the form 

cd_^ fpq^ pq^ q^ 
(7^/2 \ a ' b ' c ^ d 
The dual problem 

D* [a, b, c, d; p; q)g{z] v; a, b; p) = X(v)D* (a, 5, c', d'\p] q)g{z; v; a, b; p) (8.6) 
has a solution 

g{z; v; a, b; p) ^ T ^az^\ 6z^\ ^j^v^'^ z^^;p, q^ , (8.7) 

which is also defined up to the multiplication by an arbitrary function ^{z), (p{qz) = 
(p{z). Now it is not difficult to see that the scalar product of functions (|8.4p and 
.711 leads to the ^/-functions: 



{g{z;v;a,b;e)J{z;w;c, d\e)) 

,r ( u -iPq pq fed [cd -I \ „^ 

= V\a,b,^—v,^—v — ,— , W— W— w ■,p,q\. (8.8) 
\ y ab y ab c d \ e Ve I 

Thus, the elliptic analogue of the Gauss hypergeometric function appears to be 
directly related to the generalized eigenvalue problem for a linear combination of 
the Sklyanin algebra generators. It is necessary to note that our scalar product 
(|8.5p is different from the Sklyanin invariant measure [82]. In [80], using the latter 
measure Rosengren has built an integral representation for the elliptic Gj-symbols 
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and proved the Sklyanin conjecture on the reproducing kernel for representations 
in the space of theta functions. 

The non-uniqueness in the choices of functions / and g can be fixed by the re- 
quirement that these functions satisfy simuhaneously the equations obtained from 
(|8.3p and (j8.6|) by the permutation of p and q (because the equations >p{qz) — 
(p{pz) = ip{z) lead to (p{z) = const). This means that we introduce into consid- 
eration a second copy of the Sklyanin algebra, obtained from the first one by the 
permutation of r and 27]: 

For these two algebras the following cross-commutation relations are valid: 

SaSb = SbSa, a, 6 €{0,3} or a, fee {1,2}, 

SaSb = -SbSa, a e {0,3}, foe {1,2} or a e {1, 2}, 6 e {0, 3}. 

One can substitute in equations (|8.3[) and (|8.6p the parameterization z = e^^™/'^^ ^ 
27] ^ UJ1/UJ2, T — uj^/uj2, and to build their solutions well defined for |q| = 1 with the 
help of the modified elliptic gamma function G{u;lu). In this case the uniqueness 
of solutions can be reached by the requirement that they satisfy simultaneously to 
equations obtained from the original ones by the permutation of wi and aj2- This 
leads to another copy of the Sklyanin algebra, which is obtained from the first one 
by the transformations rj — > l/(4?7),u u/(2ry),r T/(2r]): 

'~>a — ' ' 



2') 



9 



o ; - -a+1 1 



27] 



In this case some of the generators Sa and Sa anticommute with each other as well. 

The described direct products of the Sklyanin algebra pairs can be considered as 
elliptic analogues of the Faddecv modular double Uq{sl2) ® Uq-i{sl2) [IL 84 . The 
latter double can be obtained from the second case in the limit Im(T) — > -l-oo (in 
the first case this limit is not defined) [79] , 

9. Partial fraction decompositions and determinants 

Expansions of different rational functions defined as ratios of two polynomials 
into partial fractions are used in the proofs of many identities for plain and q- 
hypergeometric series and integrals. At the elliptic level these rational functions 
are replaced by ratios of products of theta functions, and one searches for their 
expansions into sums of ratios of theta functions with the minimal number of poles. 
If the partial fraction expansion for an arbitrary rational function is a standard 
procedure, it is not so in the theta functions case. The first known relation of such 
a type follows from an identity given in [85j as an exercise. 

Theorem 7. Let 2n variables ai, . . . , a„, 5i, . . . , 5„ e C* satisfy the constraint 
rife^i — 11^=1 and Uj/ak 7^ p*^, fc e Z, for j ^ k. Then the following relation 
for theta functions is true: 



^ Y{j=i,^k^i"-k/aj-p) 
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Proof. We replace in (|9.ip n hy n + 1, denote an+i = t and substitute bn+i = 
ai . . . Gnt/bi . . . bn- After taking out of the sum the (n + l)-st term, this relation 
can be rewritten in the form [86] 

yr e{t/bk;p) _y e{tai ■ ■ ■ a,J aM ■ ■ ■ b^, p) n;^i%rV6j;p) 
^^^e{t/ak;p) ^^0{t/ar,ai---an/bi---bn;p)Yl"^j^^^^9{ar/aj;p) 

and interpreted as a partial fraction expansion over theta functions. Then the proof 
of this identity is rather elementary. For n = 2 it is reduced to the addition formula 
(jA.Sp . By induction it follows that the left-hand side can be decomposed into the 
sum 

9{tai ■ ■ ■ Un/arbi ■■ - bn^p) 



E 



Oit/ar;p) 

where the coefhcients Cr are easily found after the multiplication by 0{t/ar',p) and 
the choice t ~ Ur- □ 

Theorem 8. I87j Let 2n variables ai, . . . , a„, 6i, . . . , 6„ € C* satisfy the relations 
ttjUk, aj /a^ ^ p'', fc S Z, for j ^ k. Then the following identity for theta functions 
is true 

Proof. After the replacements n — s- n + 1 and a„^i t in (|9.3p and singling the 
(n+ l)-st term out of the sum, we obtain the partial fraction expansion of the form 



n;=i^(tof ;P) ti (^it4'-^P)UU,^kO{akaf;py 
which is easily proved by induction. □ 

These expansions into "simple" fractions for theta functions were used in the 
papers [33l [36l [86] for the proof of some exact summation and integration formulae 
for elliptic hypergeometric functions. Let us describe also another expansion which 
was used recently in |59j : 

where nLi 11^=1 = P""^- 

An elliptic analogue of the Cauchy determinant has the form 

det I — ; 5-^1 I (9.4) 

l<^.J<n \ar^e{a,zf;p) J 

±1. „\ TT .,-1;]/., .,±1. 



= ( ^y^(n-l)/2 lll<^<J<na^ '^K^f i P) ni<.<j<„ '9{z,zf;p) ^ 

Ui<Lj<n a.^^0ia^zf^;p) 
The Frobenius determinant has the form 

/ e{ta,b,-p) \ ^ e{tY['^^^a,h-p)Y{i<,<j<^ajb,e{a,/aj,h/hj\p) 
i<^j<n\0{t,a^bj;p)J 0{t;p) Ul<^,J<n0ia^bJ;p) 
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In [53], Warnaar suggested a new determinant for theta functions 

l<i,j<n \ 9{bXi,bc/Xi)n^j 

l<2<j<n i—1 

where 9{a)n — 11^=0^ 9{aq^;p). For p 0, it reduces to the Krattenthaler determi- 
nant [gS]. 

Formulae ^^-^^ are used in [IH [31 [H [HI [SHI H] and some other papers 
as auxiliary tools for proving necessary elliptic hypergeometric identities. Partial 
fraction decompositions and determinants are somewhat equivalent to each other. 
For instance, if one expands determinant (|9.5p along the last row and evaluates 
each term by the same formula ()9.5p in the smaller dimension n — 1, then there 
appears an identity equivalent to (|9.2p . Therefore formula (|9.5p follows from (|9.2p 
by induction on n, and vice versa. In a similar way, formula (j9.4p is equivalent to 
()9.3p (9T] . Applications of determinants at the level of g-hypergeometric functions 
are described, for example, in [20]. In the paper (9^, Rosengren and Schlosser 
have systematically considered determinants of theta functions on root systems (in 
particular, this paper contains a detailed list of references on this subject) and 
constructed a number of new exactly computable cases, which we skip for brevity. 
For applications of elliptic determinants to some problems of the number theory, 
combinatorics, and statistical mechanics, see [93 ] I94 j . 




10. The ELLIPTIC BETA INTEGRALS ON ROOT SYSTEMS 

10.1. Integrals for the root system C„. 

There are two different generalizations of the elliptic beta integral (|3.ip to mul- 
tiple integrals for the root system C„ (or BCn), suggested by van Diejen and the 
author [32l[35. We describe first the multiparameter integral of type I. 



Theorem 9. Let zi,...,z„ G T and complex parameters ti, . . . ,t2n+i and p,q 
satisfy the constraints \p\, \q\, \tj\ < 1 and 11^"^ ~ PI- Then 

-r-r 1 -TT nm=l^r(^'»^f ^P'g) dzi dzn 

i</i<„ n4W^) M nzf^P.q) 

T{trats;P,q), K„ = , . 10.1 

l<m<s<2n+4 ^ ' 

Proof. We consider the function 

, n T-r2n+3 t^/, ±1 n 

o(zt-C)- TT -rr n-n^i^"^' '^'^^^ 

'^"^ " liL r(z-.f 11 Vizf^.Azt^p.q) 



Y\l<m<s<2n+Z ^{tmtslP, 

where A — Y[^=i ■ For all Zi there are poles of (|10.2p in the points 
V = {i™gV, A-i9'^+V+i}„^i,...,2„+3,a,b=o,i.... 



(10.2) 
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converging to zero. The coordinates of the poles going to infinity form the set V ^ . 
Then the statement of the theorem may be rewritten in the form 

where the contom' C in an arbitrary deformation of T separating V and V^^ . 
The integral kernel p{z,t;Cn) satisfies the equation analogous to (|3.4p : 

p{z, qti,t2, . . . , t2n+3; Cn) - p{z, t] Cn) 
n 

= X! ■ • -^Q'^^t, ■■■,zn,t)- gz{z,t)) , (10.4) 

where 

Dividing equation (|10.4p by p{z, t; C„), we obtain 



Both sides of this equality are invariant under the transformation zi — s- pzi and 
have equal sets of poles (singularities at the points zi = Zj^zJ^,j — 2,...,n, 
and zi = ztp*^/^, fc S Z on the right-hand side are cancelled) with their residues. 
Therefore the functions on both sides of the equality (|10.6p differ only by an additive 
constant, independent on zi. This constant equals to zero which follows from a 
trivial check of equality (|10.6p at Zi — ti. 

Integrating (|10.4p over the variables z G C", we obtain 

I{qh,t2,...,t2n+3)-Ht)^y2\ f ~l \9^{z.t) — , (10.7) 

\JC'-ix(g-iC)xC"-* Jc^J Z 

where I{t) — p{z,t;Cn)dz/z and g^^C denotes the contour C dilated with 
respect to the zero point. 

Poles of the function pO.Sp in variable Zi converge to zero along the point Zi = 
tmq'^P^, A^^q'^p^^^ and diverge to infinity at Zi = t^q~^~°'p^'', Aq^'^p^^^^ , where 
m = 1, . . . , 2n + 3, a, 6 = 0, 1, . . .. For \tm\ < I and \p\ < \A\ the region 1 < < 
\q\^^ does not contain the poles, so that we can set C — T, deform back q^^T to T 
in (|10.7p and obtain the equality I{qti,t2, ■ ■ ■ , ^271+3) = 

Repeating almost literally the procedure of analytical continuation used in the 
n = 1 case, we find that / is a constant, which depends only on p and q. Its value 
is found by considering the limits tjtj^n ^ Ij i = Ij ■ • ■ j analogous to the n ^ 1 
case, which yields the right-hand side of (|10.3p . □ 

Formula (jlO.ip was suggested and partially justified in [52], and it was proved 
completely by different methods in [331 ESI US]- In a special limit p — > 0, it is 
reduced to one of the integration formulae of Gustafson [SS] . 
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The root system C„ consists of the set of vectors from M" of the form X{Cn) ~ 
{±2ei, ibci ±6, , i < i} I ■ , where Ci is an orthonormal basis of R". Denoting 

Zi = exp(e2), we see that the denominator of integral's kernel (jlQ.ip contains a 
product over roots of Cn of the form 

71 

QeX(C„) j=l l<i<j<n 

The root system i?C„ contains additionally the vectors {±ei, . . . , ±e„}, but the 
general rules of the appearance of these vectors in integrals' kernels are not estab- 
lished yet. The Weyl group of these systems Sn x is a symmetry of the integral 
kernel. 

The C„-elliptic beta integral of type II is built with the help of formula pO.ip 
by a purely algebraic means |33] . 

Theorem 10. Let complex parameters t,tm{m, = 1, . . . ,6),p and q satisfy condi- 
tions \p\,\q\, \t\, \t,n\ < 1, andt'^'^^'^Yfrn^i^m^pq- Then 

T{tzfzt';p,q) " nLir(t„.zf ;p,g) dzi dz^ 



T{P;p,q) 
r(i;p, q) 



\[ V{t'-H^ts;p,q)\ . (10.8) 



Proof. We denote the integral on the left-hand side of (jlO.Sp as In{t, ii, . . . , ts) and 
consider the (2n — l)-tuple integral 

[ T-r 1 fr ntor(t.zf ;p,g) 

'^"'""^ A."- ril„ r(.f .-;p,,) n r(z-;,,,) 

l<fc<Tl-l 

r(«;±4«-3/2nLii.;p,'?) dwi dw^_,dz, dz^ , 

X , „ , -I ■ E , (10.9) 

fJlT{wf\wft^^^~y^l\Lits;P,q)^i Wn-l zi 

with p, q, t and tr, r = 0, . . . ,5, lying inside the unit circle such that nr=o = 
pg. Integration over the variables Wj with the help of formula (|10.ip brings the 
expression (|10.9p to the form T^^{t;p,q)r^^{t^;p,q) In{t,ti, . . . jts) (where it is as- 
sumed that te = pqt^~^"' / Y[^=itj)- Because the integrand is bounded on the 
contour of integration, we can change the order of integrations. Then the integra- 
tion over the Zfe-variables with the help of formula pO.ip converts expression (|10.9p 
to 

r"-'(<;p,g) n ^itrts;p,q)ln^iit,t^^hi,...,t^/%). 

0<r<s<5 

As a result, we obtain a recurrence relation connecting integrals of different dimen- 
sion n: 

In{tM,---M)^^^^r^ n ntrts-,P.q)In-l{t,t''h,,...,t''h^). 
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Using known initial condition at ?i = 1 (13. ip . we find (jlO.Sp by recursion. □ 

If one expresses via other parameters and removes the muhiphers pq in the 
arguments of the elhptic gamma functions with the help of the reflection formula, 
then it is easy to pass to the limit p ^ for fixed parameters. This leads to a 
multiple q-heta integral of Gustafson |96j . A number of other limiting transitions 
in parameters leads to the Selberg integral - a fundamentally important integral 
with a large number of applications in mathematical physics [97] : 



(10.10) 



n 



l<j<n l<j<k<n 

Via + {j - l)7)r(/3 + (j - l)7)r(l + j7) 



,<^.<„ r(a + /3+(n + j-2)7)r(l+7) ' 

where Re(a),Re(/3) > and Re(7) > - min(l/n, Re(a)/(n - 1), Re(/3)/(n - 1)). 

Therefore formula (jlO.Sp represents an elliptic analogue of the Selberg integral. 
It can also be interpreted as an elliptic generalization of the Macdonald-Morris 
constant term identities for the _BC„-root system f55^. The given proof is taken from 
the work [33]. It models the proof of the Selberg integral suggested by Anderson 
[55] and represents a generalization of the method used by Gustafson in [SB] for 
proving the corresponding g-beta integral. 

10.2. Integrals for the root system An- 

Different elliptic beta integrals on the root system An have been proposed in 
papers [11[ 136] . By analogy with the C„-cases integration formulae depending on 
2ri + 3 parameters will be considered as the type I integrals. Exactly computable 
integrals with a smaller number of parameters, which are derived with the help of 
type I integrals, will be classified as the type II integrals. Let us list these formulae 
omitting their derivations. 

Theorem 11. Let \q\ < 1 and 2n + 4 parameters tm, Sm, Tn = 1, . . . , n + 2, 

satisfy the constraints \tm\,\^m\ < 1 ST — pq, where S ~ Ilm'tfi '^m o'^'^ 

/" 1 dz 

< 1 n Hz-z-i z-'z-v a) n n r(s™z„t„z-i;p,g)^ 

•^T" i<j<fc<„+i U^J^j h'P^D j^lm^l ^ 

n+2 n+2 

= l[^{Ss;n\Tt;n';p,q) H rM™;^,^), (lo.u) 

m=l k,m=l 

where Z\Z2 - - ■ Zn+i = 1 and 

A _ (p;p)So(9;9)So 



(27ri)"(n + 1)!' 



In this type I integral we have a split of 2n + 4 parameters (homogeneous in 
the C„-case) with the fixed product into two groups with n + 2 elements. This 
integration formula was proposed and partially justified in [llj . various different 
complete proofs are given in [3^ H3] . In the simplest possible p ^ limit there 
appears one of the Gustafson integrals [55] . 
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The root system An consists of the vectors X{An) — {ci — ej,i^ 
where Ci is an orthonormal basis of M"+^. These vectors he in the hyperplane 
orthogonal to the vector E = Y^^il^i s-i- Setting Zi = exp(ei — E/{n + 1)), we 
obtain zi ■ ■ ■ Zn+i — 1- The denominator of the kernel of integral (110. contains 
a product of the form 

n r(e";p,g)= ^(zt/zj, zj / z^■,p,q). 

aeX(A„) l<i<j<ri+l 

The full integral kernel is invariant with respect to the A„-Weyl group Sn+i- 
There exists an additional independent A„-integral of type I [36! . 

Theorem 12. Let \p\,\q\ < 1 and 2n + 3 parameters tk, k — l,2,...,n, and 
Sm, m = 1, 2, . . . ,n + 3, satisfy the constraints \tk\ < 1, \sm\ < 1 and \pq\ < \Stk\, 
where S — JY^=i ■ Then the following explicit integration formula is true: 

-pr ^{Sz~^z-^;p, q) 
'ttj riLl r(tfcZ-,;p, q) n^tfi r(s„2:~^p, q) dz 

""/i m=instkz-';p,q) 

n n+3 p/, N 

= 1111^^1^^^ n nSsT\-^\p.qUlOA2) 

where z\ - ■ ■ z„+i = 1. 



Here we have a split of 27t,+3 independent variables into two homogeneous groups 
with n and n + 3 elements. In the limit p ^ there appears a g-beta integral which 
was not considered in the literature until the derivation of the elliptic well 
as its plain hypergeometric degeneration appearing in the g ^ 1 limit. 

There are several A„-elliptic beta integrals of type II [TT] . We denote 



l"{t,s;p,q;An)^K^ f TT 



T{tZiZ j;p, q) 

dzi dz„ 



n+1 n+1 4 

X l[l[r{tkzT^;p,q)l[r{s,z,;p,q)^ ■ ■ (10.13) 

1 T 1 - 1 Zl Zji 

] = l fc=l 4=1 

where js^j < 1, fc = 1, . . . , n + 1, i = 1, . . . , 4, and 

n+1 4 n+1 

f^-'WtkWs.^pq, X{z,^l 

k=l i=l j=l 

Theorem 13. Under the described restrictions on the parameters the following 
An-integration formulae are true. For odd n, we have 

n+1 4 

l''{t,s;p,q;An)^ [] r(tt,tfc;p, 9) [] J] r(ifeS,;p, g) 

l<j<fc<n+l fe=l 'i=l 

T{t 2 A;p,q) i<,<„<4 
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where A = nfe=i ^fe/ /o*" e.ven n 

n+l 4 

l''{t,s;p,q;An)= r(tt,tfe;p, g) r(tfeS,;p, g) 

l<j<fe<n+l fe=l i=l 

These formulae contain only n + 5 free parameters. They can be derived as direct 
consequences of the elliptic beta integrals of type I (jlO.ip and (jlO.lip [TT] . In the 
simplest limit p ^ they are reduced to the Gustafson-Rakha (/-beta integrals 

m- 

In order to describe a different A„-integral of type II, we need ten complex 
parameters p, q, t, s, ti,t2, ts, si, S2, S3 with the constraints (fs)"~^ nfc=i ^kSk = pq 
and \p\, \q\, \t\, \s\, \tk\, \sk\ < 1- We define now the integral 

l'\t,s;p,q;Ar,)^ (10.16) 

Theorem 14. Under the indicated restrictions on the parameters the following 
integration formulae for the root system An are true. For odd n, we have 

I^\t,S]p,q]An) ^T{t^ ,s^]p,q) J]^ T{t^t4k-, SiSk]P,q) 

l<j<fe<3 

(n+l)/2 3 

X n n r((ts)^-4,sfc;p,g) (10.17) 
j=l i,k=l 

(n-l)/2 / \ 

X II iT{{tsy;p,q) Y[ T{t^-hH,tk,t^s'-'s.,Sk;p.q)\; 

j=l \ l<i<k<3 J 

for even n 

3 

l"{t,s;p,q;A„) = Y[r{t^t„s^s,;p,q) 

n/2 / 

X r(it-Hii2i3, s5-isiS2S3; 9) n r((is)^;p, q) 

j=i v 

[] r((is)^-iMfe;p,g) J] r(i^-isJi,ifc,t^sJ-is,Sfc;p,g) 



Tl/2 

(10.18) 



3 

X 

i,k=l l<i<fc<3 / 

These formulae can be derived as direct consequences of the C„-elliptic beta 
integrals of type I and II and the A„-integral of type I pO.lip [TT]. In the simplest 
possible limit p — > these integrals are reduced to the Gustafson g-hypergeometric 
integrals [M] . 

11. Some multiple series summation formulae 

There are several known summation formulae for multiple elliptic hypergeometric 
series representing multivariate extensions of the Frenkel-Turaev sum [TTl [Ml [531 
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l86l [901 flOTl [102] . A direct connection with the sums of residues of particular pole 
sequences of the kernels of elliptic beta integrals on root systems was established 
for some of them. Let us describe some of these formulae without proofs. 

Theorem 15. Let \p\ < 1 and the parameter q € <C is not equal to an inteqer power 
of p. Then the following summation formula is true 

^ i</i<„ OitMt^';p) 



j=l,...,n 



X 

.-U-1 „„-!„-! „J,-1„-1 



^ e{t]q^^^-p) e{t,tr)x, 

11 flC^2.„N 11 



9{qt-itk)N.d{qtjtk)Nk 



X 



e{qt])N, 

X 



iMn 0{qt,a-\qt,b~\qt,c-\q^+^- + -+^-~^n-\~^b~^c-^)n^ ' 
where 6 {a) X — T{aq^;p, q)/r{a;p, q) and 

q~^Y[r^^^^' ~ 1 i^^^ balancing condition), 

q^Hjtn+j = 1, J = 1, . . . ,n {the termination condition), 

with Nj = 0, 1, . . ., j = 1, . . . , n, and a = t2n+i,b = t2n+2, c = t2n+3- 

This formula was deduced in [1011 from the C„-elliptic beta integral of type I 
(which was not proven at that moment yet). Its first recursive proof was obtained 
in |86| . In the limit p — s- it degenerates to a multivariate gipY-sum, which was 
found in [inS]. 

Theorem 16. Let N ^ 0,1, . . . and the parameters p, q,t,to, . . . ,t^ £ C satisfy the 
restrictions \p\ < 1 and 

g-i^2n-2 j-jS^^ _ 2 I'^/jg balancing condition), 

q^t"^^tot4 ~ 1 (the termination condition). 

Then the following summation formula for a multiple elliptic hyper geometric series 
is true 

' e{TkTjq^^+^^ , TkT^^q^^^^^ ; p) 



0<Ai<A2<-<A„<Af l<j<k<rn\ " \' I: ' ] , ' tc ' j 

0{tTkTj)x^+X^ 0{tTkT[^)x,-X, 



TkTj,TkT^ ^■,p) 



e{qt WkTj)x^+Xj e{qt-^TkTj ^)xk-X, 



frfOirlf^A e{trT,)x, 

11 f)(T--2-„\ 11 



X 



e{rl.p) ^},e{qtv'T,)x, 



" e(qt-+=-Hi)^ ni<.<.<3 e{qt^-H-H-^)N 

f},e{qt^-^^-0Y{l^^tV^)l,Y{l^,6{qP-H^tV^)N 
Here we use the notation tj = t^P^^ , j — 1, . . . ,n. 
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This summation formula was suggested by Warnaar 53J. It follows from the d- 
elliptic beta integral of the type II [32] and was proven for the first time recursively 
in the paper f45' . 

For the root system An the type I elliptic hypergeometric sum has the following 
form: 

odr^.' M r<,M<n 

j = l,...,n 

''e{q/d,q/e)\x\ W e{t,q/h,t,qlc)x, 
_ S{q/bd,q/cd)iN\ A d{tjq,tjq/bc)Nj 

~ e{q/d,q/bcd)^N\ f \ 0{tjq/b,t,q/c)N, ' ^ ' ' 

where |A| = Ai H h A„, \N\ ^ Ni + h N,, and bcde = g^+I^L This formula 

was proven recursively in j86] . It follows also from the analysis of residues for the 
type I integral (jlO.lip \TT\ and for p — reduces to the multiple 8(/37-sum of Milne 

Omitting a number of other established summation formulae for multiple elliptic 
hypergeometric series, we describe a hypothesis from [llj. 

Conjecture. Let \p\ < 1, N — 0,1, ... , and n parameters t^. e C, k = 1, . . . ,n, 
satisfy the balancing condition Hfe^i ^fc — 9^^- Then the following summation 
formula is true: 



E 



Xi + ... + \„=N 

e{i)-N 

n even, 



0{t-'^)-N n„+l<.<,<„+3 0(i("+2)/^i.i,)-iV ' 

n odd. 



N 



It is conjectured also that this formula is related to sums of residues for the A„- 
integral (110. 13|) . For p ^ there appears a Gustafson-Rakha summation formula 
for a multiple sVT-series [lOOj . 

12. Symmetry transformations for multiple integrals 

General elliptic hypergeometric integrals of type I for the root system C„ have 
the form 

^n™''(^li ■ • ■ 7^2n+2m+4) 

n2ri+2m+4 -r^^-, ±1 n , 

^=l ^{UZj ;p,q)dzi 



11 -pi' ±1 ±1 ^ 11 T^/ ±2 \ 
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where \tj \ < 1, 



For n = we set /g™'' = 1. The integral l}^'' comcides with the ^-function, and 
I^"-' - with the eUiptic beta integral (jlO.ip . After degeneration of /n™ ''-functions to 
the level of ordinary beta integral, they reduce to the Dixon integrals ^51 1105] . In 
|34j . Rains proved the following transformation formula; 

l!rHh,...,t2n+2m+4)= J] r(M, ; p, 9) f ^ , . . . , 

l<r<.<2«+2™+4 ^ *1 t2„+2m+4 

(12.1) 

It represents a direct generalization of the third symmetry transformation for the 
function described in (|5.4p . 

In [SS], the following determinant representation has been found: 

where a^, 6j are some arbitrary parameters. For the choice = tj, bj = tn+j, 
j — 1, . . . , n, there appears the determinant of the matrix 

Tq{ti) ^Tp{tn+j) 1/^™"^" ^\qtl, ■ ■ . ,qtmptn+l, ■ ■ ■ ,pt2n;t2n+l, ■ ■ ■ ,t2n+2m+i), 

where Tq[tk) is the g-shift operator, Tq{tk)f{tk) = f{qtk)- For m = 0, we ob- 
tain thus an exactly computable determinant of univariate elliptic hypergeometric 
integrals. 

In order to prove (|12.2p , it is necessary to write 

l<*<i<«^^* l<i<j<„ l<i<j<n 

and connect both multipliers on the right-hand side with the elliptic analogue of 
the Cauchy determinant (|9.4p . Using the Heine formula 

^ / det 4'i{zj) det V'iC^i) TT d^i^i) = det [ (j)i{z)il;j{z)dfi{z) 

l<i<7i 

with (j>i{zj) = ai/9{aiZ^^:p) and ipii^j) ~ bi/9{biZ^^;q), one obtains the necessary 
result. 

The following {n + 2)-term recurrence relation takes place [33l [59] 

n+2 ^ 

X! T-rn+2 f77T±l :^9(*0^n"''(^l> ■ • ■ , i2n+2m+4) = 0, (12.3) 

where j^2n+2m+4 _ Indeed, if we write the same recurrence relation 

for the /^"'•'-integral kernel, then it reduces to identity (|9.3p after taking away a 
common multiplier. Integrating the latter relation over appropriate contours, we 
obtain ^2M- 
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The determinant representation permits one to derive another (to + 2)-term 
recurrence relation 

Enm+3<i<2n+2m+4^(*»*fe/9!P)rp /, \-l , \ n fTOA\ 

, rr flTTTT 'T-'-lW ' [tl, ■ ■ ■ ,t2n+2m+4) = ^J, (12.4) 

fc^l tkUl<^<m+2■^^k^it^/tk;p) 

where ti ■ ■ ■ t2m+2n+i — ip<lY"'^^<l- Note that transformation ()12.1|) maps equation 
(|12.3p to (|12.4|) . and vice versa. Therefore an analysis ol common solutions of these 
equations provides an alternative way of proving transformation (|12.ip [5^ . 

Let us consider symmetry transformations for other multidimensional integrals. 
For the parameters t,ti,...,ts € C satisfying the restrictions < 1 and 

^2n-2 pi^^^ tj = p'^q'^, we define the integral 

I{h,...,ts;t;p,q)^ l[ r{t,tk;p,q,t) (12.5) 

l<j<fc<8 

T{tzfzf';p,q) " ULi ntkzf;p, q) dz, 
11 r^.±i,±i-„ n\ 11 



n^r^'-^P^ 1) f=\ r(zf ;p, q) z, 



Here 



r(z;p,9,0= n (l-^iV'7')(l-^"^t^'+V+^'7'+^) 

j",fc,i=0 

is the elliptic gamma function of a higher order, which is connected to the Barnes 
multiple gamma function V,^{u\uj) and satisfies the equation 

r{tz;p,q,t) = T{z;p, q)r{z;p, q,t). 

Function (|12.5p generalizes the F-function to the type II integrals for the root 
system BCn- In [34, , Rains has proved the following symmetry transformation 
formula: 

I{ti, . . .,ts;t;p,q) = I{si, . . . ,ss;t;p,q), 

where 



Sj ^ p j = 1,2,3,4 ^ _ Itit2hti _ j pqt 



pt„ j=5,6,7,8 ' V Mil-" ~ V^sWs' |iMi.-Us,l <1. 

It describes a generalization of the key relation for the E'y-group (|5.2p . i.e. integral 
P2.5p is invariant with respect to all transformations of this group. As shown in 
|58j . the hmction I{ti, . . . ,ts',t;p,q) emerges in the quantum multiparticle model 
of [62j under certain restrictions on the parameters (the balancing condition) as a 
normalization of a special eigenfunction of the Hamiltonian. 

General elliptic hypergeometric integrals of type I on the root system A„ are 
defined as 

^n™'' (Sl, • • ■ , Sn+m+2; ^Ij ■ ■ • i tn+m+2'-: ^) 

=«n/ n ?7 — — tII n ^{sizj\tizj-p,q) 

•^T" i<,<fc<„+i r(z,^, ,z, Zk;p,q)fJl fL\ 
where \tj\, \sj\ < 1, 

n+l ?i+m+2 

j=i 
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and we set /q™'' — n;=i^ r(s;, ti;p, g). Using the result of exact computation of 
the /i°^-integral (jlO.llj) . one can derive the following recursive relation for the 
/n'"''-integrals in the variable to: 

j^m+l)j-^^^ ■ • ■ , Sn+m+3;il, • ■ • , tn+m+S', A) 

^+l•r) n\ 11 r^i-"-!/;-, , , o.s,- n o") L„ 11 



n+1 n+2 

i=i 1=1 

7fmW ^ ^ dWn 
X 7^ -'(WWI, . . . ,WW„+i,S„+3, . . . ,S„+,„+3;tl, . . . ,tn+m+2;A) 



where 



11 ~ T-rri+m+2 , T-rn+m+3 
fe=l llfc=l ^fcllz=n+3 ^1 

Q 

symmetry transformation [llj 



For TO = the /i^"* -integral on the right-hand side is computable and yields the 



InH^i, . . . , Sn+3', tl, ■ ■ ■ , tn+3; ^) = Jl ^ t„+3Sk, — ^, S„+3tfc, 

k=l V 



■p,q 



where Jlfeii ^kSk = (pg)^ and 



, n+2 

n+1 _ ''n+3 TT _ Pq 

„fj 11 ~ ^ rr"+2f ■ 

PI k=i •'^"+3llfc=i ife 



Since the left-hand side of this relation is symmetric in the parameters tk or Sfe, 
the same S'„+3 x Sn+3 symmetry is valid for the right-hand side which leads to 
additional nontrivial transformations [34 . It is convenient to denote 

n+2 

Hv) = Y[riv"^'hn+3Sk,v''+hn+3tk;p,q) 
fe=l 

X /^^^ [VSI,. . . , VSn+2,V^'^tn+3;V^^tl, ■ • • , V~^tn+2, t;"s„+3; A), 

where the arguments of In'' lie inside the unit circle, Ofc^i ^fe — Ilfc^i ■^fe — and 
V is an arbitrary free parameter (the total number of free parameters is equal to 
2n + 5). Then the derived relation can be rewritten in the form I{v) = I{v~^). 

Another type of transformations for the /^'"-'-integrals was found by Rains [33]. 
We denote T = n"=r'^^*j: 5" = lVj=r^^Sj, so that ST = and let aU 

\tk\, \sk\, \T'^^ /tk\, \S'^^/sk\ < 1- Then the following symmetry transformation 
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is true: 



n+m+2 



-^n • ■ • I in+m+2] Si, ■ • ■ , 



Sn+m+27 



X 



( 




Sri+m+2 



5^ 



) 



(12.6) 



This relation generalizes transformation (j5.3p . and there are no natural analogues 
of other iJy-reflections. We see thus that the ^-function symmetries, generated 
by different elements of the Weyl group for the exceptional root system £^7, have 
some multidimensional analogues. However, different reflections are generalized 
to different integrals whose kernels obey symmetries in the integration variables 
related to different root systems. 



Despite of a rather large volume, many problems were considered in this review 
only fragmentarily, a significant number of statements was given without proofs, 
and a number of interesting questions was not touched at all. Let us list some 
of the skipped achievements of the theory of elliptic hypergeometric functions and 
indicate several important open problems. 

Suppose that there exists a finite difference operator of the first order which 
maps given rational functions to different rational functions with a smaller number 
of poles (the "lowering" operator). In the paper [106j it was shown that this is 
possible only under the condition that the poles of these rational functions are 
parameterized by a general elliptic function of the second order, and the problem 
itself is related to the Poncelet mapping. 

We skipped description of the connection between biorthogonal rational functions 
and the Pade approximation with prescribed zeros and poles [TO t I106|, I107j . So, in 
the paper [107j it is shown that the i2Vii-elliptic hypergeometric series appears in 
the Fade interpolation tables of some functions. 

It is natural to expect that the multiple elliptic beta integrals define a mea- 
sure in biorthogonality relations for some functions of many variables, which would 
generalize the univariate relations (|7.8p . The first system of such multivariate func- 
tions, which is based on the elliptic analogue of the Selberg integral (|10.8p . was 
built by Rains [34j |35] with the help of raising and lowering operators. In certain 
limits, these functions degenerate to orthogonal polynomials of Macdonald [98] . 
Koornwinder |108j , or interpolating polynomials of Okounkov [109j (on the connec- 
tion with the latter polynomials see also [102j ). The results of the papers [34l [35] 
represent the most advanced achievements of the theory of elliptic hypergeomet- 
ric functions of many variables. Another type of generalization of the Macdonald 
polynomials to the level of theta functions was suggested in [IIQJ. 

There is a beautiful geometric interpretation of some of the elliptic hypergeo- 
metric functions in terms of the dynamics on algebraic surfaces. In jlllj . Sakai 
gave a classification of discrete Painleve equations connected with the affine Weyl 
groups. On the top of this scheme one has the elliptic Painleve equation related 
to the root system E%. In the paper [112j this equation was considered in detail 
and a reduction to the elliptic hypergeometric equation was found. Respectively, 
it was indicated that the elliptic hypergeometric series 12V11 provides a particular 
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solution of the elliptic Painleve equation. A similar role is played by the general 
solution of the elliptic hypergeometric equation [SU |58] and by the i?C„-elliptic 
hypergeometric integral of type II for some special values of the parameters |91| . 

The Pochhammer and Horn approach to the functions of hypergeometric type 
[lj[2], which we used in the case of univariate functions, is not generalized yet to the 
level of elliptic hypergeometric functions of many variables. For that it is necessary 
to learn how to solve systems of difference equations of the first order for the kernels 
of multiple series or integrals with the coefficients which are elliptic functions of all 
summation or integration variables [101 111] . This task is quite complicated, and 
all the examples, which were considered above, are built on the basis of different 
constructive ideas. On this route there appears a general problem of classification 
of all types of the elliptic beta integrals and of their multiparameter extensions to 
the higher order functions. It is expected, in particular, that there exist elliptic 
generalizations of multiple q-beta integrals for the exceptional root systems |113j . 
Let us remark also that under appropriate restrictions on the parameters all the 
integrals considered above are the integrals over polycycles. It would be interest- 
ing to consider integrals over more complicated regions of the complex integration 
variables. 

In the paper [S], a nonlinear discrete integrable system was considered and a 
self-similar reduction of the corresponding equations was suggested, which leads 
to elliptic solutions with many parameters. Using this result, a system of discrete 
biorthogonal functions was built, which are expressed in terms of the i2Vii-series. 
Biorthogonal functions described in section 7.1 represent only a particular subcase 
of these more general functions, which are expressed as linear combinations of 
several i2Vii-series and contain three additional parameters. A detailed analysis 
of these functions is not performed yet. Let us mention among open problems 
related to the general solution of the elliptic hypergeometric equation a search for 
an explicit form of the non-terminating elliptic hypergeometric continued fraction 
and a buildup of elliptic analogues of the associated Askey- Wilson polynomials 

m; 

Different generalizations of the integral transformation (j6.9p to multidimensional 
integrals on root systems were suggested in 06], but some of the corresponding 
inversion formulae are not proven yet. The problem of convergency of infinite 
elliptic hypergeometric series requires a deep analysis. It is necessary to understand 
in what sense such functions can exist. At this moment it is completely unclear 
what are the elliptic analogues of the number theoretical properties of the plain 
hypergeometric functions. It is necessary also to clarify whether it is possible to 
build nontrivial functions of the hypergeometric type for Riemann surfaces of a 
higher genus (the simplest example of such functions is given in |115j ). 

In conclusion, we can state that the main structural elements of the theory of 
plain and g-hypergeometric functions have their natural elliptic analogues. More- 
over, various "old" hypergeometric notions acquire a new meaning connected with 
the properties of the elliptic functions. Sufficiently many applications of elliptic 
hypergeometric functions in mathematical physics are known at present: in the ex- 
actly solvable models of statistical mechanics related to the elliptic solutions of the 
Yang-Baxter equation [5] and the Sklyanin algebra [551[7^[5nilll6j . in nonlinear in- 
tegrable discrete time chains , in relativistic quantum multiparticle models of the 
Calogero-Sutherland type [58j , and in nonlinear discrete equations of the Painleve 
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type |9H I112j . It is natural to expect that with time the number of such apphca- 
tions will grow and, besides that, there will appear new conceptual intersections 
with other parts of mathematics. 

I am deeply indebted to my coauthors J. F. van Diejen, E. M. Rains, S. O. 
Warnaar and A. S. Zhedanov for a fruitful collaboration and many stimulating dis- 
cussions. During the work on the theory of elliptic hypergeometric functions dis- 
cussions of various problems with H. Rosengren, S. Ruijsenaars and M. Schlosser, 
who made their own essential contribution to the development of this theory, as 
well as with G. E. Andrews, R. Askey, C. Krattenthaler, A. Levin, Yu. I. Manin, 
V. B. Priezzhev, D. Zagier, and W. Zudilin were quite useful. This review is par- 
tially based on author's habilitation thesis ^31] and lecture notes to an introductory 
course read at the Independent University (Moscow) in the fall of 2005. Some of the 
presented results were obtained during the visits to the Max Planck Institute for 
Mathematics (Bonn), to the directorate of which I am grateful for the hospitality. 
This work is supported in part by the Russian foundation for basic research, grant 
no. 08-01-00392. 

Appendix A. Elliptic functions and the Jacobi theta functions 

The periods of a periodic function are called primitive, if their linear combina- 
tions with integer coefhcients yield all periods of this function. It is well known 
that nontrivial meromorphic functions cannot have more than two primitive periods 
[117j . Functions of a real variable can have only one primitive period. These state- 
ments were used in the derivation of the expression (|2.2ip and in the construction 
of elliptic analogues of the Meijer function. The meromorphic functions f{u) with 
two primitive periods are called elliptic functions, that is there exist uji,uj2 G C, 
Im(a;i/a;2) ^ 0, and f{u + loi) = f{u + UJ2) = /(")■ 

Primitive periods of an elliptic function f{u) form parallelograms of periods and 
f{u) is determined by its values inside of them and on a pair of adjacent edges. We 
call as a fundamental domain D the interior of one of such parallelograms chosen 
in such a way that on its boundary dD there are no divisor points of f{u). Clearly, 
f(u) (as a meromorphic function) has a finite number of zeros and poles inside D. 

The integral {2ni)~^ /^^ f'{u)/f{u)du defines the number of zeros in a domain 
D for entire functions and the difference of the numbers of zeros and poles for 
meromorphic functions. Because of the periodicity, this integral is equal to zero 
for elliptic functions, that is the number of zeros equals to the number of poles in 
D. This number of zeros s (or poles) is called the order of elliptic function. The 
equality f{u) = C, where C is an arbitrarily chosen constant, is satisfied in D 
precisely s times (i.e., f{u) is a s-sheeted function). This statement follows from 
the fact that the elliptic function f{u) — C has s poles in D and, consequently, 
precisely s zeros. 

The sum of residues of the poles of f{u) in D is equal to zero, which follows from 
the equality f(u)du = 0. Therefore there are no elliptic functions of the order 
s = 1, and an elliptic function of the zeroth order is constant (Liouville's theorem). 
This gives a method of proving elliptic functions' identities: if the difference of 
functions fi{u) — f2{u) (or the ratio /i(m)//2(u)) is elliptic and contains no more 
than than one pole in the fundamental domain, then fi(u) — /2(u) = const (or 
/i(u)//2(m) = const). 
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The well known elliptic function of Weierstrass p{u\uji,uj2) jll7j has in the fun- 
damental domain one pole of the second order, i.e. s = 2. The pair [x, y) — 
(p(u), p'{u)) defines a uniformization of an elliptic curve ij^ = Ax^ — g2X — 53. El- 
liptic functions form a differential field, and any two elliptic functions f{u) and 
g{u) are related by an algebraic relation P{ f, g) — 0, where g) is a polynomial 
of its arguments. For the choice g{u) = f'{u), one sees that any elliptic function 
satisfies some nonlinear differential equation of the first order P{f, /') = 0. Tak- 
ing g{u) = f{u + y), we obtain P{f,g) = 'Ek=aPkif{u),y) f{u + yf = 0, where 
Pkifiu), y) are some polynomials in f{u) with the coefficients depending on y. Per- 
muting u and y, we see that Pkif (y) , u) = Pkif (u) , y) are symmetric polynomials of 
f{u) and f{y) with constant coefficient. The condition P{f, g) = can be rewritten 
therefore as Q{f{u), f{y), f{u + y)) = for some polynomial of its arguments Q. 
When such a condition is satisfied, one says that a function f{u) obeys an algebraic 
addition theorem. As shown by Weierstrass, a meromorphic function f{u) obeying 
such an addition theorem must be either an elliptic function or its degeneration to 
a trigonometric or a rational function. 

For a more explicit representation of elliptic functions one needs theta functions. 
Arbitrary entire functions /(u) are called (elliptic) theta functions, if 

f{u + coi) = e-^+'fiu), fiu + L02) - e™+'^/(w), (A.l) 

for some a,b,c,d€ C and Im(ti>i/cj2) 7^ 0. Replacing u by UJ2U and multiplying 
f{u) by e"" '^^'^ with some specially chosen constants a and /3, one can reach the 
equalities 

f{u+l) = fiu), f{u + T)^e^-+'f{u). 

In this appendix we use the parameterization r — uji/lu2 and denote q = e'^'^" (in 
difference from the main body of the review, where r = ^3/^2 and p = e^'^*'^). 
For a 7^ 0, the parameter b can be removed by the shift u ~* u — b/a. For a 
parallelogram D with the vertices (0, 1, 1 + t,t), we find /^^ f'{u)/f{u)du — —a. 
Therefore a = — 27ris, where the quantity s = 0, 1, . . . determines the number of 
zeros of f{u) in D. The key characteristics of a theta functions s is called its order. 

The periodicity f{u + l) = f{u) permits to expand f{u) into the Fourier series 
/(^) = X^jl-oo Cje^"^'^". Substituting it into the second equation with a — — 27ris 
and solving the emerging recurrence relation for the coefficients Cj , we find 

s-l 

f{u) = ciz' g'''=('=-l)/2(q'z^)^ z = e^'^™. 

1=0 k£Z 

The coefficients cq, . . . , Cs_i are arbitrary, i.e. theta functions of the order s form 
an s-dimensional vector space. 

If we restore arbitrary quasiperiodicity multipliers, then a theta function without 
zeros is equal to e^^*^"-', where P2{u) is a polynomial of the second order. Theta func- 
tions of the first order with one zero in the fundamental parallelogram of quasiperi- 
ods are called the Jacobi theta functions, and the functions with s > 1 are called 
theta functions of the higher level. For s = 1, it is convenient to work with four 
theta functions with characteristics 

feez 
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where the variables a and 6 take the values or 1. The standard Jacobi theta 
functions are defined as [12] : 

eMr) = Oi{u)^-ei^{u), 

02{u\t) = e2{u) e,o{u) = + 1/2), 

03(w|r) = d,{u) = eooiu) = e"-/4+™0i(u + 1/2 + t/2), 
6'4(w|t) = 6*4(71) - 0oi(w) - -ie''*^/^+"*"6'i(w + t/2). 

Note that all of them have the form X^feez''^ with h{k) = Cfc+i/c/. = q'^y for some 
constant y. Since h{k) is rational in q'', Jacobi theta functions represent a special 
class of g-hypergeometric functions. 

The 0i(u)-function is odd, 0i{—u) = —9i{u), and satisfies the quasiperiodicity 
conditions 

9^{u + 1) = -01 (m), 6i{u + t) = -e-""-2^™0i(u). 

It is related to the shortened theta function 6{z;q) — {z\q)ryo{qz^^]q)oo by the 
Jacobi triple product identity 

e^iu) = igi/Sg— «(^. q)^f?(e2-»; 5). (A.2) 

Transformation properties of the 0i-function with respect to the PSL{2, Z)-group 
of modular transformations r (ar + b)/ {ct + d), a,b,c,d£ Z, ad — be — 1, are 
determined by the relations |118| 

0i(u|t + 1) =e"/40i(w|r), Oi (^^|^) = -«V^ e""'/"0i(w|r).(A.3) 

It is convenient to use notation 9a{ui, . . . , Uk) Oaiui) ■ ■ ■ Oa{uk) and 0aix ± 
y) ■— 0aix + y, X ~ y). Then the argument duplication formula has the form 

/ 1 Y 1 + T 
6'i(2m) = —Oi [u,u+ -,u+ -,u — 

The addition theorem for theta functions, which is called sometimes a Riemann 
relation, uses products of four theta functions 

ei{u±a,v±b) - ei{u±b,v±a) ^ ei{a±b,u±v) (A.4) 



e{xw^^,yz^^;p) - e{xz^\yw^^;p) = yw-'^0{xy^\wz^'^;p). (A.5) 

The proof of this equality if elementary. The ratio of expressions standing on its 
left- and right-hand sides is a bounded function of a; £ C* (it is invariant with 
respect to the transformation x px, and it does not contain poles in the annulus 
\p\ ^ \x\ < 1). By the Liouville theorem this ratio does not depend on x, but for 
a; = w it equals to 1. 

Any theta function f{u) of the order s with the quasiperiods wi,W2 and coordi- 
nates of the zeros in the fundamental domain ai, . . . , Os can be represented in the 
form 



/(u) = e-=(-)n^iP^-), (A.6) 



fc=l 



UJ2 



UJ2 
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where P2{u) is some polynomial of u of the second order. Indeed, the function 

- afc 



9{n) = n ^ 



1 

fc=i 



UJ2 



is a theta function of the order s with the same zeros in the parallelogram of 
quasiperiods wi and U!2 as the function f{u). Therefore the ratio f(u)/g(u) is an 
entire function without zeros and poles, namely, a theta function of the zeros order, 
i.e. e^^^'^y In [118j . a self-contained theory of Jacobi forms — the functions obeying 
transformation properties similar to those of functions l|A.6p . was formulated. Note 
that all vectors fj{u),j — 1, . . . , s, of any basis of the space of theta functions of 
the order s can be represented in the indicated form with the matrix of zeros ajk 
satisfying a number of constraints (e.g., J2k=i '^jk = const). 

We call as meromorphic theta functions ratios of theta functions of an arbitrary 
finite order. It is easy to see that they define meromorphic solutions of equations 
(jA.l[) . For this we denote as ai, . . . , a„ coordinates of the zeros and as 6i, . . . , b,n 
coordinates of the poles of the corresponding function f{u) in the fundamental 
domain. Then the ratio 

is an entire function without zeros satisfying the equations /(u + uji) — e° f{u) 
and f{u + uj2) — e"^ with some a', b', c', d' , that is a theta function of the zeroth 
order e^^^"). 

Any elliptic function f{u) of the finite order s with the periods lui,lu2 can be 
represented in the form 

where C is some constant, and and bk denote coordinates of some zeros and 
poles of f{u) congruent to the zeros and poles in the fundamental parallelogram of 
periods. The following constraint should be satisfied by ak and bk'. 

ai + . . . + as = bi + . . . + bs mod lu2. (A. 8) 

This follows from the fact that both parts of equality l|A.7p are meromorphic and 
doubly periodic. Therefore their ratio defines a bounded entire function, i.e. a 
constant. The linear constraint on the values of and bk (|A.8[) . which we call the 
balancing condition, follows from the requirement of cancellation of the quasiperi- 
odicity multipliers of the 0i-functions appearing from the u —^ u + uji shift. 

Substituting expression (|X2|) in fO)) and denoting z = e^''™/'^" , tk = e^^'"''"'/"" , 
Wk = e~^'^'''"'/'^\ we obtain 



0{tkz;q) 

k=l k=l 

where the sign ambiguity appears from the factor e^™'' in (|A.2p 
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